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Abstract 

The ground-state energy of integrably-twisted theories is analyzed in finite volume. We derive 
the leading and next-to-leading order (NLO) Liischer-type corrections for large volumes of the 
vacuum energy for integrable theories with twisted boundary conditions and twisted S-matrix. 
We then derive the twisted thermodynamic Bethe ansatz (TBA) equations to describe exactly the 
ground state, from which we obtain an untwisted Y-system. The two approaches are compared by 
expanding the TBA equations to NLO, and exact agreement is found. We give explicit results for 
the 0(4) model and for the three-parameter family of 7-deformed (non-supersymmetric) planar 
AdS/CFT model, where the ground-state energy can be nontrivial and can acquire finite-size 
corrections. The NLO corrections, which correspond to double-wrapping diagrams, are explicitly 
evaluated for the latter model at six loops. 
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1 Introduction 



The AdS/CFT correspondence in the planar limit can be described by a two-dimensional integrable 
quantum field theory. The finite-volume energy levels of this integrable theory correspond on one 
side to the string energies in the curved AdS§ x S 5 background, while to the anomalous dimensions 
of gauge-invariant single-trace operators on the other side. Integrability provides tools to solve the 
finite-volume spectral problem exactly. (For recent reviews with many references, see [UEI-) 

For large volume, L, (long operators of size L), the asymptotic Bethe ansatz OH] determines the 
spectrum including all polynomial corrections in L^ 1 . In the weak-coupling limit, this result is exact 
up to L loops; but over L loops, wrapping diagrams start to contribute [5j. In the integrable quantum 
field theory, they show up as exponentially small vacuum polarization effects: virtual particles circling 
around the space-time modifies the energy levels ||6j . These effects have a systematic expansion which 
counts how many times virtual particles encircle the space-time cylinder (or diagrams wrap around). 
The leading-order (LO) Luscher correction corresponds to a single circle or wrapping. Together with 
the asymptotic Bethe ansatz, they provide an exact result up to 2L loops. The next-to-leading (NLO) 
Luscher correction corresponds to two circles and double wrapping. Including their contribution de- 
scribes the energy levels/ anomalous dimensions exactly up to 3L loops. 

For an exact description, valid for any number of loops, one has to sum up all virtual processes. For 
the ground state, this is done by the thermodynamic Bethe ansatz (TBA), which evaluates the saddle 
point of the partition function for large Euclidean times in the mirror (space-time rotated) description 
[HI EE EE El Ell E3 HI] ■ The TB A provides coupled integral equations for infinitely- many unknown 
functions, whose solutions determine the exact ground-state energy and satisfy the so-called Y-system 
relations, which is characteristic for the model and are the same for all the excited states [15]. What 
is different for the excited states is the analytical structure of these Y- functions |I6[ 1X71 ITS] . Using 
additional inputs, such as discontinuity relations [191 I2T1] and analytical structure, the Y-system can 
be turned into integral equations for excited states [HI [35], which provide the solution of the finite- 
volume spectral problem. An ultimate solution would be to replace the infinite Y-system with a finite 
T-Q system (see attempts [23j 24, 25, 26, 22j in this direction), which would lead to nonlinear integral 
equations (NLIE) for only finitely-many unknowns. 

In the present paper, we would like to analyze the ground state of the three-parameter family of 
7-deformed planar AdS/CFT theories [571 12H1 1211 1201 131]) for which we refer as 7-deformed theory 
from now on. Contrary to the undeformed or /3-deformed theories, in the most general case, no super- 
symmetry is preserved, so the ground state is indeed nontrivial and affected by wrapping corrections. 
The planar gauge theory is nevertheless ultraviolet finite and scale- invariant [32J. This is an ideal 
laboratory to test ideas directly on the ground state, which actually contains all information about 
the theory. 

The 7-deformation can be implemented in several distinct ways: in [33] it was described as an 
operatorial twisted boundary condition (the twist depends on the particle number); in |34[ [35] as a 
(c-number) twisted boundary condition and a twisted scattering matrix; finally in |36| the authors 
showed that the untwisted Y-system with twisted asymptotic conditions is consistent with the LO 
Luscher (single wrapping) correction as calculated on the gauge-theory side. In this paper, based on 
our previous work [35j , we choose twisted boundary condition and twisted S-matrix. 

We begin by analyzing in Sec [2] the effect of a twisted boundary condition on the ground state in 
general. We derive exact expressions for the LO and NLO Luscher corrections valid for any integrable 
theory with a twisted boundary condition. The LO correction contains information about the spectrum 
of the (mirror) theory, while the NLO contains the logarithmic derivative of the scattering matrix. We 
show that a Drinfcld-Reshetikhin type twist [37J of the scattering matrix does not affect the ground- 
state energy. We then demonstrate the effect of the twist in the TBA equations in general. These 
equations provides the exact description of the ground state for any finite size. By expanding the result 
for large sizes, we must recover the LO and NLO Luscher corrections. This is explicitly elaborated in 
the examples that follow. 

As a warm up in a simpler case, we analyze in Sec. [3]the 0(4) model with twisted boundary condi- 
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tions. After calculating the LO and NLO Liischer corrections, we derive the so-called raw (canonical) 
TBA equations, which contain the twist as chemical potentials. Interestingly, the twist does not show 
up in the simplified TBA equations except in the asymptotic behavior of the Y-functions. As a con- 
sequence, the Y-system is the same as the untwisted one. We solve the simplified TBA equations at 
NLO and compare with the NLO Liischer correction. We find complete agreement. 

We turn in Sec. |4]to the 7-deformed AdS/CFT model. We calculate first the LO Liischer correction. 
In calculating the NLO correction, we determine the determinant of the two-particle S-matrix S® 1 ® 2 
in all the su(2)l ® su(2)r sectors for the generic Qi and Q2 bound-state case. We then derive the 
raw TBA equations from first principles by evaluating exactly the chemical potentials originating from 
the twisted boundary condition. (For the untwisted case, the TBA equations were formulated in 
|10[ \TT\ IT2"] [T3l IT4].) The twist disappears from the simplified equations, just as it does in the 0(4) 
case. (See [38] for a general argument on this.) The twist nevertheless reappears in the asymptotic 
boundary conditions for the Y-functions. Since the simplified equations are not twisted, neither is the 
Y-system, as was anticipated by the authors of [36, 39J. Our derivation confirms their assumption. We 
then expand the TBA equations to NLO and compare with the result of the NLO Liischer correction. 
We find complete agreement again. 

We evaluate in Sec. [5] the weak-coupling expansion of the NLO Liischer correction, which cor- 
responds to double-wrapping diagrams. We explicitly compute this correction for L = 3, thereby 
obtaining the anomalous dimension of the operator TrZ 3 in the twisted gauge theory up to six loops. 

Finally, Sec. [6] contains our conclusion and outlook. 

2 Finite-size corrections of the vacuum energy 

In this section we analyze the finite-size corrections for the ground state with a twisted boundary 
condition. We consider an integrable (1 + l)-dimensional quantum field theory that possesses just 
one multiplet of particles with the same dispersion relation. The particles are labeled by a, and their 
interaction is described by the two-particle scattering matrix S^p(p±,p2), which does not admit any 
bound statesQ We are interested in the ground-state energy of a system of size L with a c-number 
twisted boundary condition in terms of the scattering data. The twisted boundary condition is defined 
by means of a conserved charge J, which commutes with the scattering matrix [J, S] — 0. The twists 
are implemented by introducing a so-called defect line on the circle. It has the effect that, whenever a 
particle of type a crosses the defect line from the left to the right, it picks up the transmission phase 
e «7 Ja ^ wnere j i s the twist angle supposed to be real. If the particle moves oppositely, then it picks up 
the inverse phase e-'T J ° . This ensures that if we formulate the Bethe-Yang equation by moving one 
particle around the circle and scattering with all the other particles and with the defect line in both 
directions, then we obtain equivalent equations. 

In deriving the finite-size energy of the vacuum with the defect line, Eq(L), we analyze the twisted 
Euclidean torus partition function from two different perspectives, see Figure [1] 

By compactifying the time-like direction with period R and taking the R — > 00 limit, the ground-state 
energy of the twisted system can be extracted from the twisted partition function as 

lim Z d (L, R) = lim Tr (e^M*) = e~ E ^ R + .... (2.1) 

R— yoo R,—>oo \ J 

In the alternative description in which the role of Euclidean time, x = —it, and space, x, are exchanged, 
the defect will be localized at a constant imaginary time t = — ix of the mirror model. It acts as an 
operator of the periodic Hilbert space of the mirror model defined by the configurations on a fixed-i 
slice. The action of this operator can be calculated from the transmission phase [5D]. In the present 
case, the operator is simply e z ' yJ , and we can evaluate the twisted partition function alternatively as 

Z d {L, R) = Tr(e-^ w V 7 ' 7 ) , (2.2) 

1 With a view to later applying this formalism to AdS/CFT, we do not assume relativistic invariance; hence, the 
two-particle S-matrix need not be a function of the difference of the particles' momenta. 
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Figure 1: Two possible locations of a defect. On the left it is located in space, and it introduces a 
twisted boundary condition. On the right it is located in (Euclidean) time, and it acts as an operator 
on the periodic Hilbert space. 



where we use a tilde ~ to help distinguish quantities in the mirror model. In the first subsection, we 
suppose that the volume L is large and expand the partition function at leading and next-to-leading 
orders. In this way, we derive the LO and NLO Luscher-type corrections for the ground state energy of 
the twisted system. Then, in the second subsection, we comment on how one can evaluate the partition 
function in the saddle-point approximation to obtain the twisted thermodynamic Bethe ansatz (TBA) 
equations. 



2.1 Large- volume expansion 

In this subsection, we evaluate the twisted partition function at LO and NLO for large volumes (i.e., 
L is large, and R — ¥ oo). This means that we keep the first two nontrivial terms in the expansion of 
the twisted partition function 

lim Tr(e-^ (fl)i e^ J ) = l + y e i7j °- g( P fc)L + V' e^-WH^HW) 1 - + . . . , (2.3) 

k,a k,l,( a <P) 

where k, I are the labels of the allowed mirror momenta p: a is the color index of the one-particle 
and (a, /3) is that of the two-particle state. The sum is taken over the distinct two-particle states. 
J is the conserved charge such that J a denotes its eigenvalue on the one particle, while J( a ,p) is its 
eigenvalue on the two-particle state. Finally, e(p) denotes the energy of the mirror particle. Clearly, 
the defect does not affect the energy levels, but nevertheless modifies the twisted partition function. 
Calculations based on the expansion of the partition function for large volumes can be found for 
boundary entropies in |41| . while for the boundary ground state energy in |42| . 



2.1.1 Leading-order calculation 

In evaluating the twisted partition function at LO, we analyze the one-particle contributions. In a 
finite but large volume, R, the momentum is quantized as 

e^ R = l -> ^-p k = keZ, (2.4) 

which is independent of the color index a = 1, . . . , N. In the R — > oo limit, the allowed momenta 
become dense, and the summation can be turned into integration. The change from the discrete label 
k to the continuous momentum variable p is dictated by the Bethe- Yang equation above as 
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Taking the logarithm of the twisted partition function, the ground-state energy can be obtained 



E${L)=- lim iT 1 log Tr{e- Ri - R)L e^ J ) 

R-^-oo 



(2.6) 



Expanding the log as log(l + x) = x + O [x 2 ) and keeping the first term, we obtain 



e*(l) = e£\l) + o( 



-2e(0)L\ 



). £«(L) = -Tr(e^) / fe~~^\ 



(2.7) 



where the color summation gives e nJa: = Tr(e l7j ), which is basically the character of the particles' 
representation. The physical meaning of this formula is clear: The finite-volume vacuum contains 
virtual particles, and they modify the vacuum energy by virtual processes. The leading volume- 
dependent process is when a particle and anti-particle pair appears from the vacuum, and then the 
particle travels around the world and annihilates with the anti-particle on the other side. Clearly, in 
so doing, it crosses the defect line and picks up the phase which, when summed up for the multiplet, 
results in the character. 



2.1.2 Next-to-leading order calculation 

At the NLO energy correction, we have to expand the logarithm of the partition function (|2.6I) to 
second order: log(l + x) = x — 2L- + (x 3 )- This will include the square of the one-particle term and 
the two-particle term. The former, however, contains a factor R 2 which would lead to a divergence in 
the R — > oo limit, and has to be canceled against a similar part of the two-particle term. We evaluate 
now the two-particle contribution and see the needed cancellation. From the remaining terms, we 
obtain the NLO energy correction. 

In calculating the two-particle term, we must first determine the allowed momenta. In very large 
volume R, the momentum quantization conditions are given by the Bethe-Yang (or, in other terminol- 
ogy, the asymptotic Bethe ansatz) equations. As the scattering mixes the color indices, we begin by 
diagonalizing the two-particle S-matrix: 

e^SZtftoPMu e iR ™e i5 »^>^ = 1 . (2.8) 

The two-particle S-matrix has N 2 eigenvalues, and we denote their phases by S^(pk,pi) for [i = 
1,...,N 2 . Unitarity implies 6/j,(pkiPi) = ~^/j,(PhPk) mod 2n. We assume that the particles are 
fermionic: S(p,p) = —I, thus 5 fJi (p,p) = n. Taking the logarithm of the equations (|2.8p for a given 
eigenvalue, we arrive at the Bethe-Yang equations 

k, 

I. (2.9) 

The fermionic nature of the particles excludes k = I; and in summing up over two-particle states, 
k > I is understood. In changing to momentum integration, it is better to reorganize the sum as 
£ fc>I /(M) = ^Ek.ifikJ) ™ |EJ(M), since the summand f(k,l) = e n^(?W+*))i is sym . 
metric. The diagonal part, — i ^ fc f(k, k), has the one-particle quantization rule (12.41) : thus, changing 
to integration as in (|2.5p the contribution to the energy turns out to be: 

E^ l \L) = \^ J f J d ±e~ 2 ^\ (2.10) 

where we used that £\ p) e J T J <°.^) = e i7j " = Tr(e^ J ) 2 . 



W~Pl- W- S fi(Pk,Pl) = 
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We now transform i J^ fe j f(k, I) into a double integral. To this end, we compute the Jacobian for 
the change of variables (fc, I) — » {pkiPl)' 



dk 


dk 


dp k 


dpi 


ai 


dl 


<>Pk 


dp, 



(2^ 



R + Sp.,k <^,i 



(27T) 



[i? 2 + i?(<W - <W] 



(2.11) 



where 5„ — dp k 5 fl (pk 1 pi) and S^j = £ M (p^ , ) . As already mentioned, the terms which contribute 
to the ground-state energy have to be proportional to R. Indeed the dangerous R 2 term 



^^iyjy f dpi I" ^P2 e -(e(pi)+e(p2))i 



2tt 



2tt 



(2.12) 



will cancel against the — ^- term of the expansion of the logarithm of the one-particle contribution. 
The second term of the Jacobi determinant (|2.11[) is proportional to the volume R, and contributes to 
the ground-state energy as 



,12 \l) = - I ^l e -^) L 

27T 



2tt 



(2.13) 



where we have used that ^(pi,^) is antisymmetric in its arguments; and that, as the twist commutes 
with the scattering matrix \e blJ , S] = 0, both can be diagonalized in the same basis. We note that 



J2e^dMp l7 p 2 ) = -z^ 1 Tr(e^ J log[5(p 1 ,p 2 )]). 



(2.14) 



In particular, this implies that if the S-matrix is twisted (a la Drinfeld-Reshetikhin [37]) with another 
conserved charge S — FSF, such that [e nJ ,F] = 0, then the finite-size correction is the same as in 
the undeformed case: 



%Tr(e^ J log(S)) 



dpi ei7J ° Tr a log(F a S a F a ) = d Pl Y, &llJa logdet(F Q S Q F Q ) 

a a 

d Pi Yl e ' nJa lo § det Sa = 9 Pi Tr ( e * 7J l °s(s)) , 



(2.15) 



where we have denoted by F a (S a ) the matrix F (S) in the subspace where J has eigenvalue J Q , 
respectively; and we have used the fact that det(F a S a F a ) = detS Q detF 2 . 

We conclude that the LO and NLO corrections to the finite-volume vacuum energy in the twisted 
theory come only from the twisted boundary condition, and are given by 



E d {L) = 4 1 \l) + 4 3 ' 1 \l)+E^(L) 

dp 



= -Th(e* J ) / ^e^^Tr^) 2 / ^ 



2tt 



dp _ 



2e(p)L 



2tt 



dp 
~2~7 



1 e -HPi)L 



dp2 
2tt 



e-^ L id Pl T, 2 {e l ~< J \og[S{p u p 2 )]). 



(2.16) 



where the omitted terms are of order of 0( e - 3c ~(°) L ), and Tr 4 () for i = 1,2 means that the trace is 
taken over the one- or two-particle states, respectively. This derivation is an alternative formulation 
of the virial expansion of the partition function in statistical physics. (See also the result for the 0(n) 
case |43j.) 

Our result (|2.16p can also be used to make the connection between the scattering description and 
other descriptions of the theory. Indeed, given an integral equation for the ground-state energy, we 
can extract from it the S-matrix by expanding for large volume to NLO. 
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2.2 Twisted TBA 



We have so far supposed that the physical volume L is large, and we have calculated the LO and NLO 
energy corrections. If the volume is not large and we are interested in the exact description of the 
vacuum, we have to evaluate the contributions of multiparticle states. This, in the untwisted case, is 
done by the TBA; and we shall now see how the derivations are modified in the presence of the twist. 

The first step in calculating the partition function is the determination of the momentum quan- 
tization of multiparticle states. This is done by solving the Bethe-Yang equations by means of the 
asymptotic Bethe ansatz (BA). Here, in addition to the physical momentum-carrying particles, one has 
to introduce so-called magnonic particles that take care of the non-diagonal nature of the scattering. 
They are useful objects, since in terms of them the scattering can be regarded as diagonal. One then 
analyzes the various "diagonal" scattering matrices and looks for bound states: i.e., complex string-like 
solutions of the asymptotic BA equations. The scattering matrices of the bound states are determined 
from the scattering matrices of their constituents. Let us label the particles (momentum-carrying, 
magnonic and their bound states) by a multilabel n; and their scattering matrices by S nm (ui, w™), 
where it™* is some generalized rapidity of a particle of type rij. Greek indices such as a will denote 
magnons only. The asymptotic BA equations for large particle numbers (thermodynamic limit) takes 
the generic form 

-l = e i ^V R lH[S nm (ut,ur), (2.17) 

m / 

where the mirror momentum vanishes for magnons p a {u a ) — 0, and 5„ n (u™,u™) = — 1. We note 
that not only the momentum, but also the energy vanishes for magnons, e a = 0. Thus, the magnonic 
equations can be inverted, without changing their physical meaning. We have to choose such equations 
which give rise to positive particle densities in the thermodynamic limit. In this limit, the partition 
function is dominated by finite-density configurations. The density of the particles (holes) of type n can 
be introduced as p n = , (p n = ^sf), where AN n (AN n ) denotes the number of particles (holes) 
in the interval (p,p + AJ5), respectively. In terms of these densities, the energy of the configuration is 

E[p] = Rj2 [ dp Pn (p)i n (p)=Rj2 [ du p n {u) e n (u) , (2.18) 

n J n 

while the entropy is 

S[p,p] = R,y^ / du[(p n + p n )log(pn + p n ) - Pn log p n - p n lag p n ] . (2-19) 
n •* 

The particle and the hole densities are not independent, and the derivative of the logarithm of the 
asymptotic BA (|2.17p connects them as 

i - _ r . ,v- 



p n + Pn - x-d u p n = j du ^ K nm (u, u')p m {u') =: K nm * p m , (2.20) 

~ m 

where K nm (u, u') = ^jd u log S'„ m (it, u 1 ). If we had inverted any of the asymptotic BA equations, then 
we would have obtained the sign-changed kernel here. By choosing the proper signs of the kernels 
for the magnons, we can ensure the positivity of all the densities. If we had started instead with the 
Drinfeld-Reshetikhin- twisted S-matrix, then S nm in (|2 . 1 T|) would be replaced by S nm , which differs 
from S nrn by constant phases; and these phases would disappear from the kernel K nm . Consequently, 
the TBA equations are independent of twists of the S-matrix, as is the Liischer correction (|2.15|) . 

We have seen that the twist does not change the energy levels of the periodic mirror system, but 
nevertheless modifies the partition function. Since the twist commutes with the scattering matrix, 
the particles of the asymptotic BA equations which diagonalize the multiparticle scatterings will have 
diagonal twist eigenvalues, too. Let us denote the eigenvalue of ijj on a particle with label n by p n . 
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The total contribution of the twist on the multiparticle state is 

n[p] = du p n (u) fin . (2-21) 

n ■* 

In terms of these quantities, the partition function can be written as 

Z d (L, R) = Tr{e-" (R)L e l ~< J ) = f JJ d[p n , p n ] e S[p,p]+^p]-LE[p] 

J n 

Evaluating the integrals in the saddle-point approximation, the minimizing condition for the pseudo- 
energies e„ = log ^ turns out to be 

e„ + fi n = t n L - log(l + e~ £m ) * K mn . (2.23) 

Once we have calculated the pseudo-energies, the ground-state energy can be extracted from the 
saddle-point value as 

E«(L) =-W pd a p n log(l + e--) . (2.24) 

Clearly the only difference compared with the untwisted case is the appearance in the TBA equations 
(|2.23j) of the chemical potential p n , which is proportional to the charge of the particle. (TBA equations 
with chemical potentials have been studied previously; see e.g. |44|.) 

As the determination of the magnons and their charges is model dependent, we work out the details 
in the following for the 0(4) model, and then for twisted planar AdS/CFT. 



3 Case study: 0(4) model 



In this section, as a warm-up, we elaborate explicitly the simpler case of the twisted 0(4) model, also 
known as the su(2) principal chiral model. We calculate the LO and NLO Liischer corrections, derive 
the twisted TBA equations, and compare the two approaches by expanding the TBA equations up to 
second order. 

The 0(4) model is a relativistic theory containing one multiplet of particles with mass m. The 
dispersion relation E(p) — y/ m 2 + p 2 can be parameterized in terms of the rapidity as 

E{9) = mcoshTrfl, p(9) = msmhn9 . (3.1) 

The particles transform under the bifundamental representation of su(2). The two-particle S-matrix 
is the simplest su(2) ® su(2) symmetric, unitary and crossing-invariant scattering matrix [431 H5] 



(0 

where 9 = 6\ — 62, and the scalar factor 



S(0) = , a .\ 2 S(6)®S(0), S(fi) =91- if, (3.2) 



So (9) = 1 2 21 [2 ' (3.3) 

yj r(i + f)r(-f) v ; 

does not have any poles in the physical strip, showing the absence of physical bound states. 

We analyze this theory on a circle of size L with a twisted boundary condition. We twist the theory 
with independent twist angles 7^ for the left and right su(2) factors, respectively: 

e * 7 j = e i 7 _J ®i+i7+i®^o = e »7- Jo ^~/+Jo = diag ^ ; q-ij g, diag ( g; ) (3 4) 

where Jo has eigenvalues ±1 on the two components of the doublet, and q — e n - , q — e J7 + . We could 
also twist the S-matrix, i.e. change S — > FSF, but this would have no effect on the ground-state 
energy, as explained in (|2.15[) . 
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3.1 Luscher corrections 

We now proceed to evaluate the Luscher correction for the vacuum (|2.16l) . As the theory is relativis- 
tically invariant, the mirror dispersion relation is e(p) = \J m 2 + p 2 , which we parameterize in terms 
of the rapidity as above: p(9) = msinh7r#. In this parameterization, the leading-order result for the 
ground-state energy is 

E^\L) =-[2],[2fem J y cosh7r0e' micosh7re , (3.5) 

where we used that 

Tr(e^ J ) = TV(e*T- J °)Tr(e^+ J °) = (q + q^Kq + <? _1 ) = [2], [2]$ ■ (3.6) 
It is useful to introduce the g-numbers 

n _ -n 

[n] q = ± q — = q^ 1 + q n - 3 + ■■■ + q 3 - n + q 1 ^ , (3.7) 

q-q 

for which [n] q — > n in the untwisted limit q — > 1. 

In the second-order correction, we have the term without the scattering matrix 

E^HL) = \WMrn J ^coshTrfe- 2 ^™-** . (3.8) 
In the other term, we have to diagonalize the two-particle S-matrix 



S(0) = S%(6)Sx(6) ® S 2 (d) , Si(fi) = S 2 (9) = -^—.S^) 



/ 1 \ 
^ 



i—i a— i 



^ /- 

V o o 01/ 



(3.9) 



The twist matrix acts on the two-particle states as 

e nJ = e ; 7 _J ^i+Jo = X ® A = diag(<f, 1, 1, <T 2 ) ® diag((? 2 , 1, 1, g~ 2 ) , (3.10) 

and commutes with the scattering matrix. The twist and the S-matrix can be diagonalized in the same 
basis, where the S-matrix eigenvalues take the form 

S = Si A ® A = Si diag(l, 1, ttl, 1) ® diag(l, 1, t±l, 1) , . (3.11) 

a — % u — % 

For the Luscher correction, we need to calculate Tr(e 47j (— idg) log S). As the scattering matrix has 
the specific tensor product structure (|3.9j) . we can write 



Tr(e l7j logS) = Tr((i ® A) (2 log S I ® I + log Si ® I + I <g> log 5 2 )) (3.12) 
= Tr(i)Tr(A)21og5 + Tr(A)Tr(ilogS*i) + Tr(i)Tr(A log 5 2 ) 
= 2Tr(i)Tr(A) log S + Tr(A) £ A t log A, + Tr(i) £ A, log A, . 

i i 

In Fourier space, the logarithmic derivatives take a particularly simple form: 

1 2t 
K o(e) = —de\ogSl(9) -> X 00 ( w - 



-l 



XW= 2^ ael ° g ^ ^ K{u) = -t\ (3.13) 
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where we have indicated the Fourier transform by tilde, and t — e 2 . The integrand of the second 
order Liischer correction is finally 

±-Tr(^ J (-id e )logS) = [2§2]^ 00 + ([2] 2 + [2]?) K. (3.14) 

In terms of these quantities, the second part of the Liischer correction is 

E { ^ 2) (L) = -[2} 2 q [2] 2 j J d6 ie - mLcosh7:ei J d6 2 coshir6 2 e- mLcosh * 92 

x {K 00 (e 1 - 9 2 ) + ([2]" 2 + [2]7 2 )A^ - 9 2 )} . (3.15) 

3.2 Twisted TBA 

Following the general procedure outlined in section 12.21 in order to formulate the twisted TBA equa- 
tions, we need to classify the particles: momentum-carrying, magnons and their bound states. We 
also have to calculate their scattering matrices; and, additionally to the untwisted case, we also must 
identify the twist charge on all the excitations. 



3.2.1 Raw twisted TBA 

In order to derive the mirror nested asymptotic BA equations, we start with an A-particle state 
consisting of down-spin particles only. We label these particles by 0. They scatter on each other as 

Soo(0) = S (e) 2 , (3.16) 

and they have the dispersion relation ea(p) = e(p). As the Jo eigenvalue of the lower component is —1 
on both su(2) sides, the chemical potential is /j,q = — 17- — 17+. We can now introduce up-spins in the 
sea of down-spins. These are the magnons, which do not change the energy and momentum, rather 
describe the polarization degrees of freedom. We label them by 1 for the right su(2) factor, and by —1 
for the left su(2) factor. Let us first focus on the positive (right) part, and denote magnon rapidities 
by u. The magnons scatter on the massive particles and on themselves as 

S m {9-u) = d ~ U+ \ , S 11 {u-u')= U - U '-\ , (3.17) 

respectively. The magnons do not have any energy and momentum ei(tt) = pi(u) = 0, but they do 
have chemical potential. Since a magnon swaps a spin from down to up, it changes the charge by 
2: /ii — 2^7 + . This means that a state with m up-spins and N — m down-spins, which contains 
type-0 particles and m type-1 particles, has Jo charge —N + 2m. Inspecting the magnon scattering 
matrices, we can conclude that a magnon and a massive particle cannot form bound states. In contrast, 
magnons among themselves can bound. Bound states in the thermodynamic limit consist of strings of 
any length M e N: 

Uj =u + i + l~ 2j , j = l,...,M. (3.18) 

We label this string as M. Clearly, the M = 1 string is the magnon itself. The scattering of the 
M-string and the massive particle can be calculated from the bootstrap, 

S OM (0-u) = f[S ol (0- Uj )= ~" + f!! . (3.19) 

3 = 1 2 
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As Sqm(0 — u)Smo(u — S) = 1) we conclude that Som(u) = Smo(u). Similarly, the magnon-magnon 
scatterings are given by 

M W 

s M m> (u - «') = n n Si1 ^ - u r ) ( 3 - 2 °) 



u' - ||M - M'|\ fu-u'- |(|M - M'| + 2)\ 2 



I u-n' + i|M-M'| / I u-u' + §(|M-M'| + 2) / 

/u - m' - |(M + M' - 2) \ 2 /u - u' - |(M + M') \ 
X "' yu-u' + i(M + M'-2) J ^u-u' + |(M + M') J ' 

These bound states have no energy and momentum cm{u) = Pm(u) = 0, while their chemical potential 
is the sum of their constituents': /j,m = 2Alij + . 

Similar considerations apply to the left excitations, which are denoted by —M. They scatter only on 
themselves and on the massive particle, such that the scattering is independent of the sign of M. The 
only difference is in the chemical potential, as the twists are different on the two sides: /i_m = 2Mry_. 

Summarizing, we have particles for any M € Z. The only massive excitation that has nontrivial 
energy and momentum has the label 0; all others are magnons. The scattering kernels in Fourier space 
have the form 

Ko°= n + l-iy K 0n =K n0 = -t n , K nm = -X_(r +m - fl"- m l) - 8 nm , (3.21) 
where t = e~^~ and n, m > 0. For the other values, we have K 0n = K _ n , K na = K_ n o, K_ n _ m = 

K-n m and K— nm — — m . 

In the general procedure, one has to invert the magnonic equations before introducing the magnon 
densities. In so doing, one obtains the "raw" (canonical) twisted TBA equations 

eo+Mo = L?o-log(l + e- eo )*X o+ J! log(l + e- £M )*^Mo, (3.22) 
cm + MM = -log(l + e - eo )*^ Af + log(l + e~ £M ') *K M >M, M^Q. (3.23) 

These equations for the untwisted (/i = 0) case reduce to those in [47], although in slightly different 
convention. 

3.2.2 Universal TBA and Y-system 

Using identities among the kernels, we now bring the TBA equations (|3.22p . (|3.23p to a universal local 
form. This means that the pseudo-energies can be associated with vertices of a two-dimensional lattice, 
such that only neighboring sites couple to each other with the following universal kernel 

a I MN = 5 MN - (K + 1) M V , s(9) = 2cQ l hn6 , (3-24) 

where I MN = 8 m +i,n + 8m-i,n and (K + l)~ h } N *{K NL + 5 NL ) = 8 M l- We also have (K n i + 5 n i)*s = 
—Ko n , which can be easily seen in Fourier space where s — t+ \-i ■ 
Let us introduce the Y-functions: 

Y M = e tM , M/0. (3.25) 



10 



We take the equations (|3.23| for Ym, act with the operator Smn — s Imn = (K + l) MN from the right, 
and use the kernel identity Kqn * {K + 1) N m = — s ^M,i- Since the chemical potentials are annihilated 
by the discrete Laplacian 

fJ-M * {sImn - S MN ) = ^(mw-i + (*n+i) - W = , (3.26) 
they completely disappear from the equations, and we arrive at 

lo g y M - Wlog(l + l«0*«. M ^°- ( 3 -27) 

Finally, we take the equations for M = ±1 and convolute them with the kernel s. We combine these 
equations with the massive equation (|3.22[) . Using the magic property of the kernel Kqo = —2s-k Kq\, 
and exploiting that /j,q + + i) = 0, we obtain the equation for the massive node 

log Yb + mZ cosh 7r0 = (log(l + Yi)+log(l + y_i))*s. (3.28) 

Thus, the twists completely disappear from the "simplified" equations (|3.27p . (|3.28p . Nevertheless, 
they enter in the asymptotics of the ^-functions as 

lim — lo g y ±M = -2<7±, (3.29) 

M— >oo 1VL 

since the kernels in (|3.23|) vanish in this limit. After all, it should not come as a surprise that the 
F-system is not twisted, 

Y M Y M = (1 + Y M -i){l + Y M+1 ) , Y ± (9) = Y(9 ± % -) . (3.30) 
The ground-state energy contains the contribution of the only massive node, 

772 f 

E (L) = -— / d6icosh7r6» log(l + Y a ) . (3.31) 
3.2.3 Asymptotic expansion 

We now make a LO and NLO asymptotic expansion of the simplified TBA equations p.27|l . (|3.28p for 

L — > oo. 

At leading order, Yq is exponentially small and the other Y functions are constant. Let us expand 
the V-functions as 

Y M = y M (l+y M ) + ..., (3.32) 

and determine all functions iteratively. The Y-system at leading order will be split into two independent 
constant Y"-systems. The solutions with the correct initial and asymptotic behaviors will determine 
the exponentially small leading-order J^o m terms of y±\. Then, in calculating the NLO j/m functions, 
we can proceed independently for the two parts. Again, the initial condition is provided by J^q, 
which appears as a multiplicative factor; while uniqueness is provided by the vanishing asymptotics 
limAf_s. 0O um — 0. The y±i obtained in this way will determine the NLO correction yo, which is needed 
for the energy correction. 

Let us now carry out these calculations. Using the fact that s * f = |/ if / is constant, we see 
from (pT2"gj) that 

log^o = -micosh7r0 + ilog(l + yi) + ^log(l + y-i), (3.33) 
where the LO constant Y"-functions satisfy the relations 

(y M ) 2 = (i + y M -i)(i + y M +i) , m ? o , (3.34) 
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as follows from (|3.27| . The solution with the correct asymptotics (|3. 291) is H 

y M = [M] q [M + 2] q , y_ M = [M] q [M + 2} q . (3.35) 

Clearly, the twist dependence reenters through the asymptotic solution. This means that at leading 
non- vanishing order 



Y w y = v/(l + yi)(l+y-i)e- mLcosh,r9 = [2] 9 [2]^e- micosh ^ , (3.36) 

which, when substituted back into the energy formula (|3.31l) . reproduces the leading-order Luscher 
correction l|3.5|) . Actually, expanding the log in the energy formula (|3.3ip to second order log(l + 3^o) = 
yo — \y^ reproduces also E^' 1 ' in (j3.8|) . Thus, we need to expand the y-functions to NLO to obtain 

the remaining Eq in (|3.15|) . 

We see from (|3.28[) and (|3.32[) that the massive node has the NLO expansion 

Y ° = y ° (' + " (irk 1 " + TT3T7"-) ) + ■ <" 7 » 

We need to calculate y±\. We expand the TBA equations (|3. 271) . keeping only the linear terms in y. 

yt + ^fe+i i + J'fe-i / 

We solve this equation by Fourier transform 

u, f -iy fl _ [fc + l]# + 3L [fc-l] g [fc+l] g 

where we have also used the result p.35[) and the identity 1 + [k — l]q[k + l] q — [k] q . Being a second- 
order difference equation, the generic solution contains two parameters. These parameters can be fixed 
by demanding that limfc_>oo Vk — and 3^o = linifc-j-o 34j/fc- The result is 

m = ffc [2] g [fc]^ 1] +2] g ([fc + 2]q - [k]qt2) y° ' y- k = §k{q ~* i] ' (3 - 40) 

which is just the deformed version of the 0(4) solution [47j . Thus, for the needed y±x, we have 



^=(^-^3^0, y_ 1= ^- M J%. (3.41) 



Performing inverse Fourier transform, 

91 = - ( K " ~w)* y " i K " - TO?) ** ■ (3 ' 42) 

Substituting back into (I3.37[) . we obtain 

Y = [2] q [2] q e- mLcosh * e (l + s*[(K 03 - [3] 9 ^ i)[2]- 2 + (K 03 - [3} q K 01 )[2]^ * [2] q [2] q e~ mLcosh ^) . 

(3.43) 

Comparing the double-convolution term with E^' 2 ^ in (|3.15p in Fourier space, we obtain complete 
agreement. 



2 The twists ■y± have small positive imaginary parts in order to suppress large-Af magnonic contributions to the 
partition function 112.221 1. 
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4 Twisted AdS/CFT 



In this section, we apply the previous methodology to the twisted AdS/CFT model. After defining 
the model by its scattering matrix, dispersion relation and twist matrix, we derive the LO and NLO 
Luscher corrections. As the model has infinitely many massive bound states Q € N, in the NLO Luscher 
correction we have a sum of the form X^Qj q 2 =i- ^ e ^ rs ^ elaborate the summand Q\ = Q2 = 1 in 
detail, and we then treat the general case, which entails detailed knowledge of all scattering matrices 
SQ1Q2 _ w e next derive the twisted TBA equations by evaluating the charges of the magnons and their 
bound states in the thermodynamic limit of the mirror asymptotic BA. The twist, just as in the 0(4) 
model, disappears from the universal equations, which lead to the untwisted y-system. We expand 
the TBA equations to NLO and compare to the Luscher correction, and again find perfect agreement. 

The AdS/CFT integrable model has an su(2|2) <£> su(2|2) symmetry. The elementary particle 
transforms under the bifundamental representation of su(2|2). For one copy of su(2|2), Latin indices 
a = 1,2 label the bosonic, while Greek indices a — 3,4 label the fermionic components of the four- 
dimensional representation. We will introduce twist in the bosonic subspace by the generator Lq, 
which has nonvanishing diagonal matrix elements: (£0)} = 1 and (£0)2 = — !• 

The symmetry completely determines the left/right scattering matrix, which has the nonvanishing 
amplitudes 




s aa = S* + S a a p - a 3 - -1 , S^p - S a p - a 4 , (4.2) 
S$ = ~\e ab e^a 5 , S% = -\e af) e ab a 6 , (4.3) 



naa c<QQ> rraa rraa (a a\ 

b aa — a 7 , O aa = a 8 , J aa = a 9 i ^aa = a 10 , (,4.4J 

where a, b <E {1, 2} with a ^ b; a, /3 <E {3, 4} with a ^ /3; and the various coefficients can be extracted 
from [48jd For Qi = Q2 = 1 we shall need explicitly only a\, since - as a consequence of some 
identities among the various coefficients - we shall be able to express the Luscher corrections purely 
in terms of it. The scattering matrix depends independently on the momenta of the particles p\ and 
P2 via 

— = e*P , x+ + -i- - x- - — = - , (4.5) 

x x + x g 

where g = \f\j (2tt) and A = g^ M iV is the 't Hooft coupling. The full scattering matrix has the form 
Sn{pi,p 2 ) = Sll {2) {p u p 2 ) Sll {m {xf,xt)® s ll(2\2){ x i> x t) • (4-6) 

sl(2) 



where SjLlpi,^) is the scalar factor 



1(2) \ 



u — u 



with a being the dressing factor. We remark that S\\ denotes actually the inverse of the AFZ S'-matrix 
[48], since we are using the relativistic convention 1 = e ipL fj • S(p,pj) , as in Section [21 instead of 
e lpL =U 3 S( P , Pj ). 

The dispersion relation can be easily expressed in terms of x^ as 

l 9 ( + 1 - , 1 



E = ~i r-F-" + Fl- (4i 



3 Indeed, 01, ... , aio are given by the coefficients of the ten terms in Eq. (8.7) in [48], respectively. 



13 



In analogy with the 0(4) model, we introduce different twists for the two su(2\2) factors, which we 
label by a — ±, 



j-yJ = e ii-L eJ7+ Lo = diag(g, q~\ 1, 1) <E> diag(g, q~\ 1, 1) , (4.9) 

where again q = e n+ , q = e ll ~ ; and "f± are related to the deformation parameters 7$ used in [29 , 30J 
by 7± = (73 ± 72)f- 

The scattering matrix has poles, which signal the existence of bound states. These states transform 
under the 4Q-dimensional totally symmetric representation of su(2|2) for any Q £ N. The dispersion 
relation of the bound states can be obtained from (|4.8[) by changing the shortening condition to 

.+ + 4-*--- = — ■ (4-10) 
x + x g 

The matrix part of the scattering matrix can be fixed |49| from the Yangian symmetry [50] , while the 
scalar factor can be determined |i5T| from the bootstrap principle. 

The mirror model has the analytically-continued scattering matrix: x^(p) — » x^(p), where p = 
—iE. Since the physical domains of p and p are different, the bound states are different, too. The mirror 
bound states transform under the 4Q-dimensional totally antisymmetric representation of su(2|2), and 
the twist charge acts as 

e^+ i0 =diag(I Q _ 1 ,I Q+1 , (Z I Q , (Z - 1 lQ). (4.11) 

The scattering matrix of the antisymmetric bound states are related to those of the symmetric ones by 
changing the labels 1 f> 3, 2 fl 4 and simultaneously nipping x ± -f-> x T inside the matrix part. Com- 
bining this with the previously mentioned notational differences, we can use the following scattering 
matrices to calculate the Liischer correction: 



$ = S?$ a (S%gfa ® S ?u\%)) > ( 4 - 12 ) 



where 



«a) = II II 5 " K > u l ) . < = U n + (On + 1 - 2j n ) - . (4.13) 

J'i=l.j'a=l 9 

and S^S?^) denotes the symmetric-symmetric bound state scattering matrix in the conventions of |49| . 
4.1 Liischer corrections 

The derivation of Section [5] is not general enough to describe the AdS/CFT problem. We have to 
incorporate two new features: the existence of fermions, and of multiple species of particles that are 
labeled by the charge Q. The fermionic nature can be taken into account by changing the trace to the 
supertrace. This is equivalent to imposing antiperiodic boundary conditions on the fermions, which 
can be implemented by an e mF twist, where F is the fermion number operator: 

STrQ^' 7 ) = Tr Q ((-l) F e^ J ) = Tr Q (e i(7rF+ ^ J) ) = STr Q (e^-' L °)STrQ(e^+- £ '°) = ([2],-2)([2]^-2)Q 2 . 

(4.14) 

Clearly, the supertrace vanishes in the untwisted q — > 1 limit. The generalization of the derivation of 
Section [5] will contain the scattering matrices S® 1 ® 2 . They arise from two-particle states with charges 
Qi and Q2- As the species are different, we should not constrain the summation on the quantization 
numbers X)fc<;i an d must keep all ^2 k ; , as they label distinct two-particle states. One can verify that 
the dangerous R 2 terms from the determinant cancel against the cross terms coming from the square 
of the one-particle contribution. Otherwise the derivation goes along the same lines as before. As a 
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final result, we obtain the LO and NLO Liischer correction as follows: 



E 



(i) 



]TSTr Q (< 



,vyu 



7^ f ^£ 
2tt 



E, 



(24) 



E ( 2 ,2) = ^ /^ e -%W /^ e -^»)iflfcSlSr 0l<fc (e^lDg^*(fc,ft))^4.17) 
Qi,Q 2 = i J *" 7F 

cf. Eqs. (|2.7p , (|2.10p and (|2.13p , (|2.14p . respectively. Here and below it is understood that = pQ i . 

In evaluating these expressions, we note that the mirror dispersion relation is defined via p = —iE 
and e = — ip. This dispersion relation can be then encoded into 



2e Q (p)L 



#2 _ 



(4.15) 
(4.16) 



rp~\~ ft nn f' 



1 



(4.18) 



where the shortening condition (|4.10p is satisfied. 

The leading Liischer correction for the vacuum (|4. 1 5[) receives contributions from each particle 



4 1) =-([2] 9 -2)([2] 4 -2)^Q 2 /f 



-e Q (p)L 



2tt 



The simple part of the NLO correction f|4. 16[) is also straightforward to compute 



E^ = i([2], - 2) 2 ([2], - 2) 2 £ Q' I ^- 2 ' eQ(m ' 



dp 



(4.19) 



(4.20) 



(2 2) 

In order to calculate the E$ ' -part of the NLO correction, we need the supertrace of the logarithmic 

(4.21) 



derivative of the mirror 5-matrix (|4~T2l) : 

STr (- e i7-io 0e i7+io log £) 

We now diagonalize the twist matrix and the scattering matrix on the same basis, 

A ® A = diag(ii, . . . , A n ) ® diag(Ai, . . . } A n ) , 
S = A® A = diag(Ai,...,A„)®diag(A 1 ,...,A„), (4.22) 



where Aj are the eigenvalues of S^^fe) > anc ^ n = 16QiQ2- Calculation similar to the one in the 0(4) 
model gives 



STr(e^ J log S) = STr (A ® A(log S^g 2 1 I + log fift* J 2) ® I + I ® log Sfjgfo ) J (4.23) 
= STr(i) STr (A) log 5^g 2 + £(-l) F * (sTr(A)i, + STr(i)A 4 ) log A, . 

i 

Using the derivative of this expression, we can express the NLO Liischer correction (|4. in the 
following form: 



E, 



.'2,2) 



J2 

Qi,Q2=l 



dpi 



2?r 



i e -« Ql (px) / dPl e -Le Q2 (p 2 ) i9 x 



2vr 



QiQ 2 (2- [2] 9 ) 2 (2- [2] 4 ) 2 log^ ( 1 g 2 (p 1 ,p 2 ) 
+ £(-l) F « (A(2 - [2],) 2 + A % (2 - [2],) 2 ) log A? lQ2 (p 1; p 2 ) 



(4.24) 
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4.1.1 NLO Liischer correction: the case Q\ = Q2 = 1 

To warm up, let us evaluate the NLO Liischer correction for the simplest Q\ = Q2 = 1 case. We focus 
on the matrix part in (|4.24[) . Performing the calculation explicitly, we obtain 

(2 - [2}g) 2 id Pl {([3], - 1) logai + log [ ai a 3 3 ((ai + 2a 2 )(a 3 + 2a 4 ) - 4a 7 a 8 )] 

-2[2} q log(a 5 a 6 - a 10 a 9 )} + (q -H- q) , (4.25) 

where (— 1) F = (1,1,-1,-1). Using the explicit expressions for the coefficients found in [48] , we 
observe the following identities 

a 5 a 6 - a w a 9 = a x , (ai + 2a 2 )(a 3 + 2a 4 ) - 4a 7 a 8 = -a\ . (4.26) 

Substituting these identities into (|4.25l) , we obtain a very simple expression for the matrix part of the 
NLO Liischer correction for Q\ = Q2 = 1 in terms of only a\, 

( 2-[2 W 2 ( 2-[2] 4 ) 2 (J?k_ + J2k_) / ^ e -^) J ^-"^^logaiCfc.fe). (4.27) 

4.1.2 NLO Liischer correction: the general case (Qi,Q2) 



Although the above approach can also be used for the cases (Qi,Q2) = (1,2), (2, 2) for which the 
explicit S-matrices are available [52], it is impractical for higher-dimensional cases. Clearly, a more 
powerful approach is needed to treat the general case. Observe from (|4.24p that the NLO Liischer 
correction involves the quantity l) Fi ^i l°g A^ 1 ^ 2 , and a similar quantity with Aj replaced by Ai. 

We exploit the fact that the sw(2|2) part of the 7+ twist e' nJ = I <£> e n+L ° involves nontrivially only 
the su(2)r factor in su(2)l <8> su(2)r C su(2|2), as is evident from (|4.1ip . Since su(2)l <8> su{2)r is the 
symmetry of the scattering matrix, we can perform an expansion in the left (sl) and right (sr) spins: 

^(-l^AaogAf 1 ' 02 = STr[(I®e^+ L «)log^ lQ2 ( SL , SK )] 

i (sl,sr) 

= E (-l) 2sR (2sL + l)[2 Sfl + l] 9 logdet^ lQ2 ( SL , Sfl ), (4.28) 

(sl,sr) 

where S'^ 1< ^ 2 (sl, sr) is the 2-particle S-matrix in the sector with left and right su(2) spins sl and sj^|, 
and we calculated the traces as Tr at (I) = 2s L + l and STr SH (e i7+i °) = (-l) 2sR {2sR + l} q , respectively. 
The sum is over all the possible values of sl and sr for the given values of Q\ and Q2- 

The problem of computing the NLO Liischer correction therefore reduces to the determination 
of det S^ 1 ® 2 (sl, sr) for general values of Q\ and Q2, and for all the possible corresponding values 
of Sl and sr. This is a formidable technical challenge, since individual S-matrix elements - and 
particularly the eigenvalues - are not known explicitly enough in general, and those that are known 
explicitly enough |52l [55] generally have very complicated expressions. Nevertheless, it turns out that 
- remarkably - these determinants have simple compact expressions, which are constructed from a 
small number of elementary building blocks. 

We propose that, with both particles in symmetric representations and Q\ , Q2 > 1, the deter- 
minants are given by the expressions in Table [T] In order to save writing, we have introduced the 
following notation 

!/o = ^. U i = \[^> U > = \1%> Uz= X { X l~\ , S Q = Ul ~ U2 ~f , (4.29) 
xj - x 2 y x 1 y xj x x x 2 - 1 Ui - 1t 2 + -f 



4 In other words, det S(sl, sr) = [], ^ifsLiSs), where Ai(sL, sr) are the eigenvalues of the 2-particle S-matrix 
corresponding to eigenstates which are also su(2)l g) su(2)r highest-weight states with given values of and sr. For 
further details, see appendix fX] As usual, the spins sl,sr are non- negative integers or half-odd integers. 
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det S^ 2 {s L ,s R ) 



?7ot r iC^2<S , 2s L +2<S , 2s L +4 • ■ • Sq 1+ q 2 -2 



||Ql2U(|Qi2| + 2),... l(Q 1 + Q 2 -2) 



(UoUtUtfS. 



HQ1 + Q2-I) 



04 

2s L + l°2s L +3 ■ 

/ r : .r ; V v -j 

Ui 



■ <->r 



1+Q2-2 



02 

i2|+2°|Qi 2 |+4 ' 



02 

' D Qi+Q 2 -2 



i(|Qi 2 | + l) J ...,|(Qi + Q2-3) 
KIQ12I - 1) > 



2s L 

(C/ C/i(7 2 )^ 2sL 5 2 4 Si+2 _ 

{U0U-JJ2) S 2sl S 2sl+ 2S2 Sl+ 4 : - 



i+C 



(C/oC/lt/2) 4 (^) 21 ^Q 12 | + 2^|Q 12 | + 4- 



' °Qi+Q 2 -2 



U UiU 2 5'|Qi 2 |+2S'|Q 12 | + 4 

{u Q u 1 u 2 f s 2 2 slsl . . 



.5, 



Qi+Q 2 -2 



UQ1 + Q2) 
|(Ql + 02-2) 

|(Qi + Q 2 -4) 

KIQ12I + 2),... 
IIQ12I 96 
KlQial - 2) > 



i(Qi + Q 2 -6) 



o3 

• °Qi+Q 2 -2 



= Q 
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Table 1: The values of det S® 1 ® 2 (sl, sr) for Qi,<5 2 > 1 and for all possible su and Sz, where 
$21 = [q = ±1) and = Qi — Qj- See (|4.29|) for further notations. 



where a;* are the parameters of the bound-state representation, and Uj ± = ar~ + Note 

that there are only three possible values of the right-spin, namely sr — 0, 5, 1, as 2s r counts the 
number of fermions in the basis of the Hilbert space. If at least one of either Qi or Q 2 is 1, then the 
corresponding results are collected in Table [5J A brief account of how these results were obtained is 
presented in Appendix [A} 

Substituting the results from Tables [1] and [2] into (|4.28[) , and carefully simplifying the resulting 
expression, we obtain 

^(-l) F '^logA? lQ2 = [3] q {QiQ 2 logU U 1 U 2 +K Q ^) (4.30) 

i 

-[2],((4QiQ a - Qi - Q 2 )log[/ + 2Q 2 (2Q 1 - l)logE7i 

+2Q 1 {2Q 2 - 1) log U 2 + (Q a - Qi) log f7 3 + 4/C QlQ2 ) 

+[l],((5QiQa - 2Qi - 2Q 2 )logC/ + Q 2 (5Q! - 4) logtfi 

+0i(5Q 2 - 4) log CT a + 2(Q 2 - Q x ) logC7 3 + 5/C QlQ2 ) , 

where we have defined 

Qi-l Qi-j-i 
IC QlQ2 = («2 - Qi + 2j + 1) £ log^Q 2 -Q 1 +2 J +2fe ■ (4.31) 

J=0 fe=l 

In deriving the result (I4.30P , we have made use of the fact that log Sq is an antisymmetric function 
of Q (i.e., log^-Q = — logiSg, up to an irrelevant additive constant), which in particular implies that 
JQQ1Q2 — JQQ2Q1 ^ w e emphasize that (|4.30[) holds for any Qi,Q 2 € N. An analogous result can be 
derived for ^2A—l) Fi A4 log A^ 1 ^ 2 by replacing q — > q in (|4.30p . 



17 



sr 




det SQ^{s L ,s R ) 



U UiU 2 



HQi + Q 2 -2) 



|(<3i + 
KIQ12I 



Q2-I) 
- 1) > 



KQ1 + Q2) 
§|Qi 2 |^o 

KIQ12I-2) >o 

= Qi2 (Qi = Q 2 



1) 



Table 2: The values of det S Qi ® 2 (sl 7 sr) for all possible s_r and sl if either Qi or Q 2 is 1. 



Thus, in order to calculate the Luscher correction, we have to plug (14.30)) into the formula (14.241) : 



(2.2) = j2 Q1Q2 r —< '■ f°° —<■ 



1=1 



2tt J_ x 2tt 

x id fl { (2 - [2]^) 2 [3], (Q1Q2 log E/0E/1E/2 + /C QlC 



-[2], ((4QiQ 2 - Qi - Q 2 ) logf/o + 2Q 2 (2Q! - 1) logf/i 

+2g!(2Q 2 - 1) log U 2 + (Q 2 - Qi) log U 3 + 4/C QlQ2 ) 

+[1], ((5QiQ 2 - 2Qi - 2Q 2 ) log U + Q 2 {$Qi - 4) log U x 
+Qi(5Q 2 - 4) \ogU 2 + 2(Q 2 - QO log[/ 3 + 5/C QlQ 

-QiQ 2 (2 - [2],) 2 (2 - [2] 4 ) 2 log5 s ^ 2 (p 1 ,p 2 ) i . 



(4.32) 



si (2) 

We shall compare this result to the TB A output in Section 14.4.41 
4.2 Twisted TBA equations 

In [191 |2"UI [T5] , the authors derived the TBA equations for the AdS /CFT model with the most general 
chemical potentials. Hence, the TBA equations for the 7-deformed theories correspond to some special 
cases. However, since we must determine precisely the charges/chemical potentials of the various 
excitations in terms of the deformation parameters, we now briefly sketch the derivation. 

In order to derive the TBA equations, we have to recall the various types of excitations (both 
massive and magnonic) and their scattering matrices; and we must calculate their twist charges. We 
label the fundamental massive particle as Q = 1, corresponding to the (33) label of the fundamental 
representation. The S'-matrix of this kind of particles is in fact given by (|4.7[) and they can form 
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bound states for any Q with string- like complex roots denned like in (|4.13j) . We label such a massive 
composite particle by Q and the scattering matrix of such particles is ([4.131) . Since the twist charge 
acts trivially in the (3,4) subspace, the massive particles are not charged: /j,q = 0. 

We now focus on the magnonic excitations. They encode the color su(2|2) structure of the scat- 
tering, and come in independent left and right copies. We first consider the right part. We label a 
magnon, which introduces label 2 in the sea of massive 3-particles, by y. It scatters trivially on itself, 
but nontrivially on the massive particles 



Si v (u,y) = S Qy (u,y) = Y[Sx y (u j ,y). (4.33) 



Q 

n 



The twist charge of the y particles is fi y — — ?7+. 

We can also introduce the label 1 in the sea of 2-particles. These particles are labeled by w. They 
scatter nontrivially only on the y particles and on themselves: 

Syw(y, w) = S-i (v(y) - w) , S ww (w,w') = S 2 (w - w') , (4.34) 

where v(y) = y + y , and S n (u) is defined as in (|4.29|) . namely 

S ^) = — t- (4-35) 

" + a 

The twist charge of these particles is [i w — 2ij + . 

As the scattering matrix S yw (y, w) has a difference form in the variable v{y) = y + y , we might 
use the parameter v instead of y. The inverse of the relation, however, is not unique. Defining 
y-(v) = i(v — iV4 — v 2 ) with the branch cuts running from ±oo to ±2, we can describe any y with 
3m(f/) < for v <E [—2,2]. Clearly y+(v) = y-(w)" 1 describes the other ?sm(y) > case; and in the 
scattering matrices Si y which depend on y, and not on v, we have to specify which root is taken. As 
a consequence, we have two types of y particles y\S with S = ±; and the scattering matrices split as 
S\ y {u,y) — ¥ Si y \s(u,v) := Si y (u,y$(v)). Clearly, the y\S magnons scatter on the momentum bound 
states as S Qy \ 5 (u,v) = JJj S ly \ S (u j: v). 

Let us now focus on the magnonic bound states. Detailed investigation showed [53J that v and w 
particles can form bound states for any positive integer M. It consist of 2M w-particles f±(M+2-2j) = 
v ± (M + 2 — 2j) ~ for j = 1 , . . . , M with yj = (yZj ) * , and M w-particles with synchronized parameters 
Wj — v + (M + 1 — 2j)jj f° r j — 1) • • • ) M. The scattering matrix of the v\M particle with all other 
particles is simply the product of the scatterings of each of its individual constituents 

M+l M Af-1 

S v \Mi(v,q)= ]J S y\-*( v M+2-2j,q) Y[ s wi(wj,q) S y \ + l (v M -2j, q) ■ (4.36) 

The twist charge of the bound state simply sums up to (i v \m = 2M(— + M2vy+ = 0. 

The w-type particles can form bound states among themselves: an A-string of w-particles can be 
formed as Wj = w + (A + 1 — 2j)|. The scattering of the A-string with any other particle is 

JV 

Sw\Ni(w,q) = Y[ S wi (wj,q) , (4.37) 
i=i 

while the twist charge is fi w \N — 2Nij + . 

We summarize the various scattering matrices and chemical potentials in Table [3] 
Once we know the chemical potentials, we can calculate the kernels 

Kjj, (u, u') = ^-8 U log Sjj. (u, u') , (4.38) 

Z7TI 
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v\M 2 


w\N 2 


y\$2 




Qi 


SqiQi 


Sq 1 v \m 2 


1 


SQi y\s 2 





v\M x 


Sv\Mi Q 2 


S v \Mi v\M 2 


1 


c 

^'v \i\ y\o 2 





w\N x 


1 


1 


S w \Niw\Na 


S W \N ± y\S 2 




v\5i 




I Si v\ M 2 


Sy\8i w\N 2 


1 


— ?7 



Table 3: Scattering matrices of the various particles and their chemical potentials for any of the two 
su(2|2) wings. 



and write the TBA equations one by one. To ensure positive particle densities, we have to invert the 
equations for v\M and for y\—. The equation for the massive nodes then read as 

£q 2 = Ll Q . 2 -\og{l + e-^)*K QlQ2 +Y, teg{l+e-<\M)*K v \ MQ2 -8\og{l + e- e y\*)*Ky\5Q 2 ■ (4.39) 

ot=± 

Note that for particles of type v\M and y\6, we must include contributions of the two sw(2|2) copies, 
which we denote by a = ±. The remaining equations are valid for the two su(2|2) factors separately, 
so we omit the a index: 

e v \M = - log(l + e~ £Q * ) * K Q2 v \ m + log(l + e" £ "iM' ) * ^ |M , y{M _ $ i og (i + e -%|s ) * K y \ 5v \M , 
e w \ N = -fi w]N -log{l + e~^^ N ')*K wlN , w \ N -Slog(l + e~ e ^ s )*K v \ SwlN , (4.40) 
e y \ S = i'tt- /ij,| a - log(l + e" e «2) * A" Q2 y |,5 + log(l + e " e »i M ) * K v \ My \ s - log(l + e~ e *°l w ) ★if t( ,|iv 1 /|« ■ 

Once these equations are solved, the ground-state energy can be obtained as 

£ (Z) = - £ / + e ~ £Q2 ) • ( 441 ) 

In [TO1 [20] the authors analyzed the TBA equations with generic chemical potentials, and formu- 
lated the requirement under which the Y-system remains unchanged. Our chemical potentials, which 
correspond to 7-deformations, satisfy their requirement. 

4.3 Universal TBA equations and Y-system 

The TBA equations can usually be brought into a local form. As already remarked, this means that 
the pseudo-energies can be drawn on a two-dimensional lattice, such that only neighboring sites couple 
to each other with the universal kernel 

s I M n = S MN - (K + 1)^ N , s(u) = - — g7rM , (4.42) 

where I M n = $m+i,n + 5m-i,n ail( l (K + 1) MN * (K NL + S NL ) = 5ml- To simplify the notation, we 
introduce the following Y-function^l 

Y Q = e-^, Y v \ M = e^\", Y w{N = , Y s = -e^ . (4.43) 

Clearly, we have two copies for the magnonic F-functions: Y3 M ,Y*< N) Y" where a = ± refers to the 
two su(2|2) copies. Acting with the operator (|4.42|) on the TBA equations (|4.39|) . (|4.40|) . and using 
kernel identities such as (K + l)7f N * Kn — sSm,i &s well as the special properties of the chemical 

5 To compare with [TJUl], we note that Y w{N = Ygfa, Y V \ M = Y^f vw and K& Q * = K^§ 2 AF . Also, K n (u) = 

- , , _M 

2^-£^ logS , n (?i), where S n (u) is defined in 14.351 1; its Fourier transform is K n = sign(n)tl n l, t = e s , 
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potentials fJ, w \N—i — 2fj, w m + ^,«|Ar+i = and fi w \\ = — 2/z y , we arrive at their simplified form [14] , For 
later purposes, we write the simplified equations for v\M and w\N magnons, and a useful combination 
(hybrid) of the un-simplified equations for Q and y particles |16] 



log Y Q2 = -Li Q2 + log(l + Y Ql ) * (K^ + 2s * K%~ hQ2 



E 

a=± 



1 - 



log f 1 + Yf {1 j * s *K yQ2 + log(l + n?g a -i) * s - log j— * 



1 



i 

1 1 _ Y 7 * „ 

log 

Z 1 yo 



+ o 1 °S( 1 -^)( 1 



log YfY£ 
log — 



1 + 



It! 1 



logKj 



|M 



2 ov y aM y 

- log(l + Y Q2 ) * Xq 2 + 2 log(l + Fq 2 ) * 1 * s + 2 log I 
log(l + r Q2 )*X Q2y 



i - r 

log(l + Ym+x) * s + I M n log(l + FJijv) * s + £mi log —*s . 



-k S . 



1 -Yl 



l°g 1T|m = Imn log(l + y^|jv) * s + 5mi log — j-*s , 



(4.44) 

(4.45) 
(4.46) 
(4.47) 

(4.48) 



where in the convolution ★ we integrate over the interval [—2, 2] only. To conform with part of the 



literature, we have renamed some kernels KMP = K, 



K% M = K f 



\M 



K 



yQ 



K 



v\-Q 



K, 



y\ + Q> 



Kq v — Kq y \_ — Kq y \ + . The ground-state energy is given by summing the contributions of the massive 
nodes only: 



E (L) 



E 

Q=X 



J ^-d u p Q log(l + Y Q ). 



(4.49) 



Evidently, as in the case of the 0(4) model, the chemical potentials and so the twists completely 
disappear from the simplified equations: They show up only in the asymptotics of the Y W \ N functions, 
as limjv-i.00 logl^i^v = — MiuliV- It follows that the F-system relations are not modified by the twists, 
as was supposed in [35]. Equations (j4.48[) - (|4.49[) together with the asymptotic prescription give the 
complete solution for the finite-size energy of the twisted AdS/CFT model for any coupling g. We now 
check this solution against LO and NLO Liischer corrections. 



4.4 Asymptotic expansion 

We now expand the simplified TBA equations to leading and next-to-leading order. We expand any 
K-functions as 

Y = y(l + y + . . . ) . (4.50) 

We solve iteratively these equations similarly to the 0(4) case: At leading order, all the massive 
nodes Yq are exponentially small, which splits the Y-system into two independent subsystems which 
have constant asymptotic solutions. These constant values then determine the LO exponentially small 
expressions for J^q. At NLO, we obtain linear integral equations for the y corrections of the two 
subsystems, whose initial values are provided by the asymptotic J^q functions. The solution of the 
linearized equations determine the NLO correction for the massive nodes yQ, which provides the NLO 
energy correction. 



4.4.1 Leading-order expansion 

At LO, the massive J^q functions are exponentially small, and we can neglect the convolutions involving 
all log(l + Yq). The magnonic y±, y„i M , 3^ijv functions are constants. From (|4.46p . we see that 
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ya = ya It tnen f n ows f rom ((07)| and (|4~48|) that the equations for y^ M and y%, N are the same 
as those for one of the wings of the 0(4) model ()3.27j) . From the asymptotic behavior, we see that the 
solution for v\M is the same as in the undeformed model, while the solution for w\N is that of the 
deformed model: 

y% M =M(M + 2), y+ lN ^[N] q [N + 2} q: y- [N =[NUN + 2} i . (4.51) 
Since 1 ★ s = | , the equations (|4.45l) for y± can be solved as 



where we have further streamlined the notation by defining 

[n] + = [n] g , [n]- = [n] q . (4.53) 

The sign in (|4.52|) can be fixed by the last equation in (|4.40|) . and is consistent with the vanishing of 
the ground-state energy (|4.56[) in the undeformed (q,q — > 1) limit. We now use that HK v q = 1 (see 
(6.12) in [13]) to write 



log^Q = —L e Q + \ £ [log(l + X Q |i) + log(l + ^ Q _x) + log(l - 3^) (l - 3^) 



(4.54) 



Using the asymptotic solution (|4.51[) , (|4.52p . we obtain the leading-order result for Yq 

y Q = (2- [2} q )(2 - [2} q )Q 2 e~ L ^ . (4.55) 
Substituting back into the energy formula (|4.49D , the LO correction reads as 

E (L) ~ e£\l) = J2 / T y Q = "< 2 ~ ®M 2 ~ ^ E <? / Y e ~ L ~ €Qm > ( 456 ) 

Q=1 J ^ Q=l J 77 

which agrees with the result (|4.19[) that we obtained from the Luscher calculation. 
4.4.2 NLO expansion 

Expanding the energy formula (|4.49[) to NLO, we obtain 

*M = -T.f£ 108(1 + YQ) - - E / g W + VQ) + E / 1^ . (4-57) 



Q=l 



The quadratic term nicely reproduces our previous result (|4.20p for Eq , since using again (|4.55p 
gives 

4 2A \L) = ^ffhyl = h2- [2] q f{2 - [2] q f Q 4 f fe~^) . (4 . 5 8) 

Q=l 77 Q=l J 7F 

In order to evaluate 

4 2 ' 2 \L) = -f: f^y Q y Q , (4.59) 
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we must first calculate tjq . This will be given by the solution of the following linearized set of TB A 
equations: d 



VQ2 



Q = ± 

k iq, , y--y+ , + 



-★if. 



yQ 2 



Vv\M 

y w \N 



2{y+ - 1) V2 2(3^ - 1 
2 (K\iVv\i ~ K\ivl\) * s - y Q2 * k Q2 + 2y Q2 * k^ 1 * s , 

* S - ^M+l * S + S M 1 ~ x + • 



A w\N-lVZ\N-l + A Z\N+lVZ\N+l) *S + 6 N l 1 _y a *S 



(4.60) 

(4.61) 
(4.62) 

(4.63) 
(4.64) 



where 



A a = 



y"\ M _ M(M + 2) 



A a — 

Ji w\N — 



w\N 



1 + v a , 



[iV]a[A r + 2] a 
[N + l]l 



(4.65) 
is only in 



We start with the equation (|4.64|) for y™i N - The difference between a = + and a 
the asymptotics (|4.5ip . (|4.52p . Since one equation can be obtained from the other by interchanging 
q <-> q, we do not write out explicitly the a index. Replacing y + — y_ in (|4.64p with the contributions 
from the massive nodes (|4.62[) . and using the explicit form of the asymptotic solution, we obtain an 
equation similar to the one for the 0(4) case: 



where 



Dw\N 



Cl 



[iV-l][iV+l] [N + l][N + 3] . 

Vw\N-l-\ [jV + 2] 2 yw\N+l) * s + °NlCl* S , 



[Nf 



[2] 



-,yQ * K Qv , K Qy (u, v) = K Qy { Ul v) {Q(v + 2) - Q(v - 2)) 



(4.66) 



(4.67) 



[2] - 2' 

and Q(v) is the standard unit step function. We solve the difference equation in Fourier space. We 
use that s = (2 cosh ^) _1 = (t + 1^ 1 )^ 1 where t = e~~ . The solution which decreases for large N (to 
respect the asymptotics of Y w \n) and is compatible with the Sn,i term is 



Vw\N 



cit f [N + l] 
[2] ^ [N] 



N-l 



[N+l] 
[N + 2] 



N+l 



(4.68) 



We now analyze the equation (|4.63[) for y^\ M - This difference equation is not the same as for the 
undeformed 0(4) model, as it has inhomogeneous terms, 



Vv\M 

where 



( {M -1){M + 1) (M + l)(M + 3) i 

-Vv\M-l H TTFi o^9 -Vv\M+l * S - y M+ l *S + 5 M \C2 * s , (4.69) 



M 2 



(M + 2) 2 
2 



C2 



[2] -2" w ~" Ww ' 
Taking the Fourier transform, we obtain the difference equation 



(4.70) 



{t + t~')y v 



(M-1)(M+1). , (Af+l)(M + 3). 



M 



M 2 



-Uv\M-l 



Vv\M+l — yM+l + SmiC2 ■ 



(M + 2) 2 

6 We note that in |17l there is an erroneous term in Eq. (2.7): — YR * s should be instead —Yq * Kq, as in 14.611 1 



(4.71) 
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The general solution with two arbitrary parameters Ai and A 2 reads as 



V ^ M - {-^T* ~ M + 2 ) 

/ M + 1 X _ M M + l i-M 
\ M M + 2 



^ y, y fc+ lt- fc - 2 (t- 2 fc-fc-2) \ 

v 1 S r 2 -l) 3 (fc + l) J 

/ A y, y fc+ it fc - 2 (t- 2 (fc + 2)-fc) ' 
V 2 £1 (t- 2 -i) 3 (fc + i) , 

The parameters can be fixed from lhxiM-s-oo y v \M — and from the M = 1 term as 



(4.72) 



Ai=t 2 (—c 2 -A 2 ) , A 2 = ]T 



y k+ it k - 2 (t- 2 (k + 2)-k) 



t (t- a - 1)° (A + 1) 



(4.73) 



The NLO hybrid equation for j/q 2 is (|4.60[) ; we plug into it the equations (|4.61|) and (|4.62[) . and 
obtain 



VQ 2 



U 1 



w\l 



-t/^i! * s*^q 2 - — — -y w]1 * s*K, 



yQ 2 



y+ 



A-v\Q-2~lVv\Q 2 -l * s 



(4.74) 



Since y v u and y w \\ can be expressed in terms of J^q, we see that the solution for uq 2 has the general 
form 

VQ, = y Ql * + y Ql *M Q i Q *. (4.75) 

Consider the first term y>Q 1 * K®^® 2 . It is easy to see that its contribution to the integrand in the 
energy formula (|4.59[) 

y Q , y Q , = y Q , (y Ql * x^a? ) , (4.76) 

with yQ 2 given by (14.55(1 . matches with the "scalar part" of the integrand of the Liischer correction 

Ep > in (1424)1 . We now proceed to analyze the remaining contribution in (|4.75p . and show that it 
gives the "matrix part" of the integrand of the Liischer correction. 



4.4.3 NLO TBA correction: the case Q x = Q 2 = 1 

To warm up, let us evaluate the NLO correction for the Q\ = Q 2 = 1 case; thus, we calculate A4 11 . 
In so doing, we can freely put 3^q 2 = for Q 2 > 1. The corresponding solutions read as 



Uw\l 
Vv\l 

y+ - y- 
y+ + y- 



[2] 



[2] 


-2 




2 


[2] 


- 2 


yi 


*K 




2 


[2] 


- 2 



yi+K ly i[K 



[3] 



(4.77) 



y x -kKty-k ( Ki - -K 3 



3 [3] 



2 [2] 



1 1 



yi*K ly i[ l--^)Ki-[-- 1 z ! )K3)*s-y 1 *K 1 + 2y 1 *K£*a 



2 [2] 



It is convenient to substitute these solutions directly into (|4.60p . i.e., 



yx 



yi*Kit {2) + 



a=± 



(4.78) 
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Using the explicit form of the asymptotic solutions, one can see that the terms involving the convolution 
with K3 completely cancel. Exploiting further that Ki y *Ki = KH (which can be shown using relations 
from Sec. 6 in [13 1), we arrive at 



[2] 



2([2] Q - 2) 

[3]a - 3 



(3^1 * K ly *K x +y x -k K 1 iK yl - 2^i * Kll -k siK yl 
2 



3>i * Kll * s * K yi 



([2]a - 2) 2 
This expression further simplifies to 

[2]e 



[2] Q - 2 



y>i * K ly is * if. 



(4.79) 



2/i 



2([2] Q -2) 



y 1 -k(K ly *K 1 + i<r 1 *ir„ 1 )+— — * s*K yl - k 1v * s * Kll) 



[2] a -T 



In the second term, using K\ v *Ki = KH and KH = K\*K y \, we can write 

Kll * siK v i - Kiyis * Kll = Ki y *(K\ *s - S* Kl)*Kyl = . 



(4.80) 



(4.81) 



as both s and K\ depend on the differences of their arguments, and therefore their convolution is 
commutative. In the previous term in (|4.80[) . we can obtain 



1 



K ly iK x + KiiKyi = ^—d Ul log 



x l — x\ x\ — l/x 2 ui — U2 — 2i/g 
Xi -l/x~2 Xi - x~2 u 1 -u 2 + 2i/g 

X I " Xr-, I X 



= -d Ul \og[ 1 2 _ \^ = ^ log 01 («!,«;,), (4.82) 



7X1 



X-\ X<-) 1/ X 



where we have used identities from Sec. 6 in (T3] and Eq. (3.7) in P3J, and we have recalled the 
definition in (|4.1[) of a\. The final expression for the Q — Q' — 1 contribution to the energy (I4.59P is 
therefore given by 



yi 



yi * Kll(2) + yi * M 11 , M u = ^-d Ul log oi( Ul> «a) 



[2]c 



27r -«i-o ,v ^~'>^2 - [21 Q ' 
which completely reproduces the result (|4.27p obtained directly from the Liischer correction. 
4.4.4 NLO TBA correction for any Qi,Q 2 

We now consider the general case. Let us recall the result (|4.74|) for j/q 2 



(4.83) 



VQ 2 



=y Ql *(K? l $ 2 +2s*K%- 1 > Q2 )+ E 



A v\l 



A' 



aQ 2 



+^"|Q2-iy"lQ2-i * s 



a=l,2 L y+ 

2(y% - 1) 



•y a _ I 1 - 4r 2(3^2 - 1) ^ 2 2(3>£ - 1) 



(4.84) 



and analyze it term by term. Since we have already checked in Section (|4.4.2ll the matching of the 
first term with the scalar part of the Liischer result, we start by considering the second term of (|4.84[) . 
which can be rewritten as 



2y Ql * a * Kg 1 - 1 '** = 2y Ql * K Q ^x * siK yQ2 + 2y Ql * E A Q2-Qi+2,+i * a , (4.85) 

3=0 
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where we used the property s*Kq = Kq*s, valid for any Q. Now we consider the terms in the square 
brackets of (|4.84|) . again suppressing the index a. Using the solution (|4.72p for M = 1 and taking its 
inverse Fourier transform, we can express the first term as 



= 23p]( 2 - [2] \T- Kl 

+ l^*[(Qi-l)KQ 1+ i-(Qi + l)K Ql -i]^*3iKyQ a , (4.86) 

where the term in the second line can be rewritten, by using the identity (K n +\ + K n —\ + nS nt ±x8)*s = 
K n) as 

2 „ Qi-l 

- 9 _ [2] * * s * K yQ2 + q 1 ( 2 - [2]) ^ Ql * KQl * KyQ2 ■ 

The if 3 contribution in the first line of (|4.86[) cancels, as we have already seen in the Qx = Qa = 1 
case, with the successive term in (|4.84|) 

- y W ^ i y w \i*siK yQ2 = - ( 2 _ 1 [ 2 ])2 y Qi *K QlV *(K 3 - [3]K!) * slK y Q 2 , (4.87) 



while the terms with if 1 give 

[3] -3 



y Ql * K Qiy iKi * siK v Q 2 . (4.88) 



(2-[2])2- 

Summing this contribution to the first two terms in the second line of (|4.84|> . we obtain 
2 



2 -[2; 



(yQi * K Qiv*s * Ki*K y Q 2 - y Ql * K QlV iKi * siK y Q 2 ) 



[21 2 

r^Qi * ^Qi-i * S *^Q 2 + o — f^T^Qi * K Qiv* s * K Q2-i , ( 4 -89) 



2- [2]" V1 V1_1 yV2 2- [2] 

where we used the identities K^® = KykKyQ + Kq-% and = KQ y iK\ + Kq-%. Q As already 
noticed for the case Qi — Q2 = 1, the first line in the expression above vanishes because K\-ks — s-kK\. 
The successive two terms in the second line of (]4.84[) give 

[ " yg 1 **a ui logo? 1<,!, (ui,«a) J (4.90) 



2^(2 - [2]) 

where we used the identity (j4.82[) generalized for any Qi, Q2, 



K Qiy iK Q2 + K Ql *K yQ2 = -d Ul log ( - 1 - J_ Q2 J j =--d Ul \oga^ Q >( Ul ,u 2 ). (4.91) 

Let us turn to the last and most complicated term. Using the inverse Fourier transform of (|4.72j) for 
M = Q2 — 1, we can write it as follows 

A v \Q 2 ~lVv\Q 2 -l * S = -r^Y * ^ k(k~ Qx) [{Q 2 + l)K Q2 -Q 1+2 k-l - (Q2 - l)KQ 2 -Q 1+ 2k+l\ * s 



yQl *KQiy 

Qa(2-[2]) 



* [(Q a - 1)^q 2 +i - (Q a + 1)^q 2 -i] * s , (4.92) 



7 The latter identity is reported in footnote 4 of |16k the former can be derived analogously using equations (6.19) 
and (6.39) in [13 . The same equations, together with (6.14), can also be used to obtain 14.851 . 
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where the second line can be expressed as 
2 



2 -[2] 



y<3i * K Qiv* s * Kqz-I 



1 



! (2-[2]) 



y<3i * K Qiy* K Q 2 



(4.93) 



Now, taking into account that in summing over a the first term in (|4.92|) gets a factor 2 and the other 
terms get similar coefficients with q — » q, we can sum all the contributions above to get the solution 
for i/q 2 for generic values of Qi, Q 2 : 



h-2 



^ * K ?ii2) 2 + 2 E K Q ^ Ql+2j+1 * . 

I 3=0 



d Ul logaf lQ2 (Mi,u 2 ) 



7H 



E 

a=± 



flu, log a? l0a (tii.ua) , d Ul loga£ 2Ql (u 2 , Ul )* 



27rig 2 (2- [2] Q ) 



27riQi(2- [2] £ 



2 Ql_1 1 
E fc ( fc ~ ^ t^ 2 + 1 )- ft: Q2^Qi+2fe-i ~ (Q2 ~ l)K Q ^ Ql+2k+1 ] *s\ 



(4.94) 



where we used the following identity 



K Qiy* K Q 



-d ui log X \ Qi ^ ^ i^i— 



27TJ 



T 9 ul loga^ lQ2 ( Ul ,u 2 ), (4.95) 



its hermitian conjugate (recall that x(u)* = l/x(u*) in the mirror kinematics) 



(4.96) 



and a£ lWa (ui ) u 2 )aS' aW1 (ii2,ui)* = o? l0a (ui,tta) 



convolutions involving the universal kernel s(u) in (|4.94|) as 



Moreover, we can write the sum of the two 



1 Q1-1 

— — V fc(fc - Qi) [(Q 2 + l)iYQ 2 _ Ql+2fe _i - (Q 2 - l)X Q2 _ Ql+2fe+ i] * s 
Q1W2 



Q1-2 



(4.97) 



where we used the definition f|4.31[) of /C^ 1 *^ 2 . Remarkably, despite the long computation, the final 
expression for yq 2 is quite simple and reads 



2/Q2 



• log a^ 2 u 2 ) + — log a^ 1 (« 2l Ul )* 



7 S ;(2) 



QiQ 



_£<3iQ2 _21ogaf lQ2 (wi,w 2 ) 



y 1 

4l(2-[2] a ) 



a=± 



1 

0^ 



1 



(4.98) 



8 Actually, identities 114.821 1. 114.911 1. and 14.961 1 are valid up to vanishing derivatives d ui log \ j — ?p, d ui log \j \q 2 and 
du! log 2 _q 2 , respectively. 



-Q2 



27 



Substituting this result, together with the result (|4.55j) for D^Qd into the formula (|4.59|) for the energy 
correction, we obtain 

x id Pl I (2 - [2] 4 ) 2 [[3], (-Q1Q2 log «? lQ2 + ^ QlQ2 

-[2], (-4QrQ 2 log af lQ2 - Q x log a® lQ2 - Q 2 \oga® 2Ql * + 4/C QlQ 

+ [1], (-5QiQ 2 log a? lQ2 - 2Qi log a^ lQ2 - 2Q 2 log a* Ql * + 5IC QlQ2 
+{q O (?) 

+QiQ 2 (2- [2] ? ) 2 (2- [2]j) 2 logS%ffi(pi,p2)\ ■ (4.99) 

Finally, through the following identifications 

a? lQa = (UoUtUz)- 1 , a^ Q2 =U UlU 3 , af Ql * = U f) U 2 l U^ 1 1 (4.100) 
we find full agreement with the result (|4.32p from the Liischer computation. 

5 Weak- coupling expansion 

In this section we calculate the weak-coupling expansion of the ground-state energy of the twisted 
AdS/CFT model. In order to perform the weak-coupling expansion, we use the parameterization 



which follows from (|4.10l) and (|4.18[) . At leading order in g , and so at weak coupling, we have 

x-=tli9. + 0(g), x + = ^JL^ + 0(g 3 ). (5.2) 
9 p + iQ 

5.1 LO contribution, single wrapping 

The LO correction can be calculated from (|4.56[) by using the expansion of the exponential term 
appearing in J^q: 

00 „2(L+j) 

e ~Le Q (p) 9 /ggx 

~h (p 2 +Q 2 ) L+r 

In particular cq = 1, while the higher-order terms can be easily generated with Mathematica. Using 
the fact that ^f^ n \z) = § ^{j^l+i , we perform the integral in (I4.56P by residues 

The summation over Q gives rise to a series of ^-functions: 

4 1) (L)^-(2-[2] 9 )(2-[2],)£c J 2 1 - 2 ^) ( 2 f +p_~ 2 ) C 2(i+ ,)-3.9 2(i+j) . (5.5) 

i=o 

This result is exact up to g 4L where the NLO Liischer correction starts to play a role. We evaluate 
the leading <? 4i -order contribution of the NLO Liischer correction in the next subsection. 
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5.2 NLO contribution, double wrapping 

The simplest term of the NLO correction comes from (|4.58[) and contains y>g . Its contribution at order 
g AL can be calculated using eq. (|5.4[) to be 

E£ A \L) = (2 - [2] 9 ) 2 (2 - [2] 4 ) 2 2-«* ( 2L 1 1 ) ' ( 5 ' 6 ) 

The most complicated term is i?g ' (L). We have to evaluate (|4.59|1 based on the solution given in 
(|4.98j) . The twist dependence comes in two distinct ways as: 



4 2 < 2) (L) = (2-[2U 2 (2-[2]^ 



A(L) + B(L) 



[2} q -2 [2] 4 -2 



g 4L - (5.7) 



We first calculate B(L) for any value of L. The weak-coupling expansion of the functions a 2 and 

Q2Q1* ■ u 
a 2 are given by 

d Pl \oga^ Q *(p u p 2 ) = 0(g 2 ), d i;i loga^(p 2 ,p 1 r = ^^ + 0(g 2 ) . (5.8) 

Pi 1 Wi 

We substitute these results into (I4.98|) and then into (|4.59p , we perform the integrals as in (|5.4I) . and 
sum up the independent terms to obtain: 

B{L) = -2 1 -" ( 2 r L " 2 )( 2 M C 2i -iC 2i -3 ■ (5.9) 



L- 1 J \ L 

This gives the complete answer for the given (2 — [2] 9 )(2 — [2]q)(4 — [2] g — [2]g) dependence of the 
double-wrapping correction at leading nonvanishing order for any L. 

We now proceed to calculate A(L). It acquires contributions from the first line of (|4.98[) . which we 
denote by A s u 2 -\, A/c and A\, respectively, 

A{L) = A sl{2) (L) + A K {L) + A X {L) , (5.10) 

where 

A sli2) (L) = £ Q\Ql f^±e- L ^ [^e- L ^\d p Jo g S%$(pi,p 2 ), (5.11) 

A K (L) = 2 Y, f ^e- L ^^ f ^e~ L ^^hd Pl lC^(p u p 2 ), (5.12) 

ML) = -2 £ Q 2 Q 2 /'|ie-^^) /|i e -^fe) l a, i loga?^(p 1 ,p 2 ). (5.13) 
In order to compute Ai, we expand a® 1 ® 2 to leading order in g: 



9^ loga? l02 (pi,p 2 ) = + 0( 3 2 ) . (5.14) 



Substituting the result back into f|5 . 13[) gives 



ML) = -2 1 "" ( )[ L ) C2L-2C2L-3 ■ (5.15) 

Observe that the transcendentality of A(L) and B(L) are different. It seems the deformation 2 — [2] 
carries transcendentality 1 . A similar effect was observed already in [5"H [55] . 
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To calculate A s [( 2 ), we have to expand the logarithm of the dressing factor logS^5 2 (pi,p 2 ) in the 
mirror-mirror kinematics. According to [T2] , it has the structure log S®lS 2 (pi , p 2 ) — ~ log ^u{2) (-Pi >P^)' 
2\ogY, QlQ2 (p 1 ,p 2 ). Hence, we can write 



^-d Pl logS«g 2 (pi,p 2 ) = -K QlQ2 - 1% IogE« l «»(pi,pa) • 
Explicitly performing the weak-coupling expansion of (6.14) in [51] . we obtain (see (|B.4j) ) 



id §x logE QlQ2 (pi,p 2 ) = -- 



VX1 - + *Q0) - ^(1 + - Pa) + Qi + Qa)) + cc 



(5.16) 



(5.17) 



where tp(x) — d x (logT(x)) is the polygamma function. The su(2) scalar factor results in 



K su(2) = K QlQl 



1 

47T 



^(*(pa - Pi) -Qi + Qa)) + ^(1 + \{i(P2 -Pi) ~Qi + Qa)) (5-18) 
-ip(\(i(p2 - Pi) + Qi + Qa)) - + ^(*(Pa - Pi) + Qi + Qa)) + cc 



Finally, 



id Pl K. QlQ2 = 



4(Qi - 1)Q 2 + ((Qi - Q2) 2 + (pi - P2) 2 ) x 



+ -(i(pa - Pi) - Qi + Qa)) - ^(^(«(p 2 - Pi) + Qi + Qa))) + cc. 



(5.19) 



(5.20) 



Denoting the contributions to A s i( 2 ) by Ay, and A su ^ 2 ^ we have that 

A sm (L) = Ax(L) + A su(2) (L) , 

where 

~ V Le «^ 2) i%log£ 0l ° 2 (pi,p 2 ), (5.21) 



A S (L) = -2^ Q?Q 



2 / d Pi e -L eQl fe) MPa 



2tt 



2tt 



A s „ (2 )(i) = ^Q^/fle-^^/"|a e -^ 2 (p 2 ) 2 ^ Qi02 ^ )#2) , (5>22) 



Using methods explained in Appendix [Bj we evaluated the integrals by residues. To demonstrate 
the structure of the result, we write out explicitly Ay, (see (|B.10[) and (|B. 13JI ) : 



^ 21 ( l L_ 1 2 )c2 i -3EE / 



L-l 



£ £ 

Qi,Q2 ii ,.72=0 



L -!- ,/ 1 — I \ 2- 2L + 2 



-2L + 1 



31 y(i-i-ii)! 

L+j 2 -n 2-^ 



y(£-i-i)i 

(-Qi) 2 " 1 '-''' 1 



(-QO^-V^-^COi + i) 



-T T (-Q2) 2 " L ^ 2 ^ (2i ~ J1 ^ 2 " 2) (Qi + Qa + 1) ■ (5.23) 



J2 J (L-l-j 2 )\ 
These terms can be expressed in terms of multiple zeta values (MZV) by rewriting 



^ n \q + 1) = (-ir +1 n!(c(n + 1) - J2 r 71 - 1 ) ■ 

3 = 1 



(5.24) 



9 For n = 0, one has to replace £(1) with je- 
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and performing the sums explicitly. We will, however, not pursue this calculation further here as we 
did not find an explicit answer for generic L. The integrals can be evaluated similarly for A su ^ 2 ) and 
Ax. with a similar structural final result, although some care must be taken to the Q\ — Q2 dependent 
term for Qi = Qi- In the next subsection, we present explicit results for the smallest nontrivial length: 
L = 3. 



5.3 Specific calculations for L = 3 



The LO wrapping correction (|5.5II for L — 2 is divergent, as we have for j = the term C2L-3 = Ci- 
Similar observations were made in |55[ 138] . We therefore focus now on the first nontrivial case, namely 
L = 3. The LO correction for this case goes as follows: 



4 1} (3) = -(2 - [2],)(2 - [2],) - ^C55 8 + ^Ci 9 U 

The simple double-wrapping contribution (|5.6[) at leading order is 

^ 2 - 1) (3) = (2-[2] 9 ) 2 (2-[2y 2 I |- 4 
(2 2) 

In calculating the term Eq : (3), we recall from (|5.7p that 



15 



945 



3465 
~256 



C99 1 



C79 1 



(2,2) /on _ 



(3) = (2-[2]^(2-[2]^ 



From (|5.9| . we have 



A(3) + S(3) 



15 



[2], - 2 [2], - 2 



^(3) = -^CsC 5 . 



We calculated the contributions to j4(3) one by one. The simplest is 

15 



A x (3) 



256 



C3C4 



(5.25) 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



as follows from (|5.15D . In the more complicated terms, we calculated the integrals by residues as 
explained in Appendix [Bj Then, in summing up the expressions, we employed the following strategies: 

• We performed the sums analytically by replacing the polygamma functions with harmonic sums 
using (|5.24p . and then rearranging all the sums into MZVs. These MZVs could then be expressed 
in terms of elementary ones, which contained only products of simple zetas with transcendentality 
less than or equal to 7. 

• Alternatively, for terms involving polygamma functions depending on Qi + Q2, we replaced the 
polygamma functions with their integral representations 



On.O- 



■l) n 



dt. 



(5.30) 



and performed the summations YIq x q 2= i explicitly. The remaining integral over t could be 
evaluated numerically with very high precision (100 digits), and the result could be expressed in 
terms of products of zeta functions (and the Euler constant je) with the help of the online MZV 
calculator, EZ-Face. 

• Finally, for polygamma functions depending on Q\ — Q2, we evaluated the sums numerically 

as q 2 =i = ^Sq! Eq2=i + ^(Qi-Q 2 )— 1 ; anc ^ a S am expressed the result in terms of zeta 
functions using EZ-Face. 



■'EZ-Face is documented in 1561 . and can be accessed at http://oldweb.cecm.sfu.ca/projects/EZFace/index.html 
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We found the following results: 



81 21 441 

As(3) = 1024 C3C4+ 512 C2C5 " 2048 C7 - 
A fo\ 9 y y 315 



9 3 91 63 

4c(3) = -— C?-^-CsC4-— C2C5 + — Ct« (5-31) 
256 S3 1024^ S 512^ ^° 512 s7 



By summing up, we obtain the total A contribution 



9 , 2 63 



A(3) = A E (3) + i4„ (a) (3) + ^c(3) + A x (3) = -— C| - . (5.32) 

Thus, the total anomalous dimension is 

E (3) = 4 1) (3)+4 2 ' 1) (3)+4 2 ' 2) (3) + ... 

, n ,„, „„ M , / 3 . « 15, o 945, in 3465, 12 \ 

- -(2 - [ 2 y (2 - [%) - C3 / - -c 55 8 + ^c 75 10 - w c, 9 12 + . . . 



-(2 - [2],)(2 - [2] 4 ) ([2], + [2], - 4)) ^C 3 C 5 5 12 + ■ • 



+(2 - [2] 9 ) 2 (2 - [2Uf (- JLtf + ^Cv) 9 12 + ■ ■ ■ , (5.33) 

where we recall that 2 — [2] 9 = 4sin(^ !l ) 2 and 2 — [2]^ = 4sin(2^-) 2 in terms of the deformation 
parameters j± = (73 ±72)5, as in our case L = 3. 

The result (|5.33|) is indeed the total anomalous dimension, since the vacuum energy does not receive 
any contributions from the asymptotic Bethe ansatz. Remarkably, even though at intermediate stages 
of the computation there appear terms involving even zeta functions and Euler's constant 7^, all such 
terms finally cancel. 



6 Conclusion 

We have computed the NLO finite- volume correction to the vacuum energy in twisted AdS/CFT 
by two apparently independent approaches: Liischer (|4.32|) and TBA (|4.99p . The fact that both 
approaches yield identical results provides a strong consistency check on the AdS/CFT S-matrices and 
TBA equations that have been developed in the literature, as well as on the final result. This result is 
expressed in terms of a double infinite sum of contributions from the infinitely-many types of massive 
mirror bound states. Our computations check the complete (both horizontal and vertical parts of the) 
F-system, and go beyond the five- loop calculations presented in |17[ \57\ [58], which checked at the 
single wrapping order only the vertical part. 

Our result is valid for any value of the coupling constant. However, by making a weak-coupling 
expansion, we have obtained a prediction (|5.33p for the anomalous dimension of the operator TrZ 3 
in the twisted gauge theory up to six loops. It should be possible to check this prediction directly in 
perturbation theory by taking into account both single-wrapping and double-wrapping diagrams. To 
our knowledge, this is the first complete computation of double wrapping in the literature. It may be 
interesting to investigate also the strong-coupling limit. 

The key results needed for the NLO Liischer computation were the determinants of the (untwisted) 
AdS/CFT S-matrices in all the su(2)l <g> su(2)r sectors, presented in Tables [1] and [5J The simplicity of 
these results suggests that they may have some group-theoretical formulation. In particular, it should 
be possible to find a general proof, presumably based on sw(2|2) Yangian symmetry. 

It would be interesting to extend our analysis of finite-size corrections in twisted AdS/CFT, which 
has so far been restricted to the ground state, to excited states beyond the LO result of [3S]. It would 
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also be interesting to understand the origin of the divergence of the LO and NLO results for L = 2, 
which was already noticed in similar contexts in |55l 138] . Finally, one can now begin to contemplate 
triple and higher wrapping. 
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A Determinants of S-matrices in the su(2)l <8> su(2)r sectors 

We describe here how we obtained the results in Tables Q] and [2] for det SW 1 ^ 2 ) (s Lj sr), the determi- 
nants of the AdS/CFT S-matrices in the su(2)l <8> su(2)r sectors, which enter into the NLO Liischer 
computation. Our straightforward approach was to explicitly compute these determinants for small 
values of Qi and Q 2 (up to 8), and then infer the general pattern. 

For the cases (Qi,Q2) = (1, 1), (1, 2), (2, 2), we used the explicit S-matrices from [52] to directly 
compute the eigenvalues. For the cases (Qi,Qa) = we used results from [59]: from Eq. (56) 

there, it follows that (up to the overall factors), 

de , S (i,« ( fcI !) = „«, 

d « s "- a, <!4> - (5 

de,s<>.« ( e^,i) . do t (;? 

detS<«)(£±±,0) = o! = l, 
detS «,«.(^i,0, . |±idet 

detSC«(^,0) = ^a\l (Ad) 
One can verify using the explicit values of a\ that 

a% = idetf a I a l )=U U 1 U 2 











«2 










o§ 


4 


4 















1 



Qdet 11) = ^r^det 4 4 4 \ =U*U X U*U; 



4 4 J Q-l 
l 



3 ■ 



a\l = U UlU 3 , (A.2) 



11 We note a couple of typos in appendix B of 1591 : a| should not have the factor ai+z+ in the denominator; and a| is 
missing an overall minus sign. 
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where the notation is defined in (|4.29|) . It is then easy to see that the expressions in Table [2] are 
consistent with the results (|A.1[) . (|A.2[) . 

For general values of (Qi,Q2), we made use of the formalism developed in [55]. As an example, 
let us consider the case (Qi,Q2) — (2,3). Since the state of a single Q-particle (the 4Q-dimensional 
totally symmetric representation of su(2\2)) has the su(2)l <S> su(2)r decomposition 

Q n Q-l l Q-2 n . 

x V° + V— x V? +V— x V°, (A.3) 
the decomposition of the corresponding 2-particle states can be obtained from the tensor product 

(y 1 xV° + V* xV* + V° x V°^j ® (yi x V° + V 1 x vi x , (A.4) 

where in this appendix we denote by x the tensor product of the su(2)l and su(2)r representations. 
For concreteness, let us focus on the computation of det S' 2,3 '(l, 5). The tensor product in (|A.4[) 
can be decomposed, by the Clebsch-Gordan theorem, into a sum of irreducible representations of 
su(2)l (£> su(2)r. In this decomposition, there appear four representations with (sl,su) — (1, |), 
which are the relevant ones for computing this determinant. These four representations come from the 
following channels: 



1 : 

2 : 

3 : 

4 : 



(y 1 x V°) <g> (V 1 x V*) 
(yh x V^j ® fv* x V°) 



(v° x v°) ® (y 1 x v* 

(yk x V*) ® (yi x V°) . (A.5) 



The corresponding highest-weight states (sl, sr)) with sl = mi, = 1 and sr = m,R = \ are 

given (up to an overall normalization factor) by 

l# 3) (l^)) <* |0,1>?-]1,0>?, 

I# 3) (i,~)) « |o,i)?, 
I# 3) (i,~)) « |i,o)?, 

|# 3) (l,i)) « |0,1>? - ]l,Q>i x , (A.6) 
respectively, where the states \k, I)] 1 are defined in [49]. It is convenient to introduce a basis \a) of 



k + l = h 




\ei) = 


|o,i>? 


N = 


|0,1)? 


|es> = 


10, l)^ 1 


|e 4 ) = 


|i.o)? 


|e 5 ) = 


|i.o>? 


lee) = 


11,0)? 



(A.7) 

Although these states are orthogonal, they are not normalized^ Indeed, defining 

m = (ei\ei), (A.8) 



We are grateful to M. de Leeuw for pointing this out to us. 
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it readily follows from the definitions of the states [49] that here = (2, 2, 1, 6, 2, 2). An orthonormal 
basis is therefore given by 

|gi) = — j=|ei) , (e 4 |ej) = <% . (A.9) 

The S-matrix acts as 

S|e i >=^|e i )C/ ji . (A.10) 

3 

Numerical values for the coefficients Uji can be computed using formulas in 49 , for given numerical 
values of momenta pi,p2, coupling constant <?, and representations Qi,Q2CJ Hence, we can obtain 
the corresponding coefficients Uji in the normalized basis 

Uji = (ej\S\ei) 

= M Uji y.u.i.iuu,/. (a.ii) 

' 1 k,l 

where we have introduced the diagonal matrix Mij = y/Tfyfiij. A useful check is that the matrix Uji 
(unlike Uji) is unitary. 

We use (|A.6[) to express the highest-weight states 1,( ^ 2 \sl, s r )) in terms of the normalized 

l4 0l ' 0a) (»L,«fl))=£|gi>Ciz, Cl/ EE(e~ l |4 Ql ' Q2) ( SL ,^)), (A.12) 

i 

where the states themselves are normalized, 

(4 QuQ2 \sL,s R )\Tpf 1 > Q2 \s L ,s R ))=5 I j. (A.13) 
We can finally construct the S- matrix in the (sl, s r ) sector, 

S<tf' Q '\8 L) 8 B ) = (^ Q2 \ SL , SR )\S\^ J Ql ' Q2 \ SL , SR )) =J2^iU ijC3J . (A.14) 

id 

Another useful check is that the matrix Sj < ^ 1 ' < ^ 2 \sl: s R ) is also unitary. Computing numerically the 
determinant of this matrix 

detS^^Hs L ,s R ) = det(s^' Q2 \s L ,s R )) , (A.15) 

we find for the case in question (namely, (Qi,^) = (2,3) and (sl,s r ) = (1, i)) that the result 
coincides with (t/ot/it^) 4 S3, in agreement with Table [TJ Other cases (Qi,Q2) and other sectors 
(sl, s R ) can be treated in a similar way. Note that sectors with s R = 1, |,0 are constructed with 
states of type I, II, III, respectively. After some effort to accumulate results for sufficiently many cases, 
the general pattern summarized in Tables Q] and [5] became evident. 

Before closing this section, it may be worthwhile to frame the problem that we have addressed 
here in a more general context. Consider an S-matrix (solution of the Yang-Baxter equation) that 
is invariant under a group G, which here is su(2)l <£> su(2) R . As is well known (see e.g. [60l l61]). 



13 We note that version 1 in the arXiv of |49l contains a number of typos, most of which are corrected in the journal. 
However, some typos remain in the latter. In particular, in (5.14): Qij = bidj — bjdi. Also, in A~ 1 in (5.17): in the (2,2) 
element of the big matrix, should be instead ; and in the (2,1) matrix element, the sign in front of [M+ (Zi — /2V2] 

should be plus instead of minus. Finally, in (A. 8), the formulas for 61, , 64 should have sign plus instead of minus; and 

the formulas for d,2 and d-j, should not have i in the denominator. We are grateful to G. Arutyunov and M. de Leeuw for 
correspondence on these points. 

14 We use the convention that the determinant of a number (i.e., a 1 X 1 matrix) is the number itself. 
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the S-matrix is a matrix S ab defined in the tensor product of two vector spaces V a (g> V b in which 
representations IT a and ITf, of G act. F^l 

[Il a {g)®Tl b (g),S ab ] = 0, geG. (A.16) 

The representation space decomposes into a sum of irreducible representations of G parameterized by 
highest weights A&, which here are (sl,sr), 

V a ®V b = J2 v ( A k)- (A.17) 

Since the S-matrix is G-invariant (|A. 16|) . it has the corresponding spectral resolution 

S ab = 5> fc P Afc , (A.18) 

k 

where Pa,, is a projector onto the irreducible subspace V(Ak). 

In the seminal work [6Q\ on the construction of rational S'-matrices, it was essential to assume that 
the Clebsch-Gordan series (IA.17|) is multiplicity free (i.e., a given irreducible representation appears at 
most once), in which case pk in (|A.18[) is a scalar. For AdS/CFT, the decomposition (|A. 1 7|) is unfor- 
tunately not multiplicity free: the Clebsch-Gordan series contains multiple irreducible representations, 
as we have seen in the example (|A.5|) . Hence, pk becomes anrxr matrix, where r is the multiplicity 
of the corresponding irreducible representation with highest weight A&. In the AdS/CFT case, pk is 
the matrix that we have defined in (|A.14| . The problem of explicitly determining this matrix can be 
quite complicated even for rational S- matrices, see e.g. [ST|- In the present work, we have restricted 
to the problem of computing its determinant. 

B Details of the weak-coupling expansion 

B.l Weak coupling expansion of the dressing phase 

The dressing phase in the mirror-mirror kinematics is given by |51j 

-ilog£ QlQ2 (2/i,2/ 2 ) = $(y+y+)-$(y+,%)-$(yf,y+) + $(yr,y 2 -) 

+\ [-*(yt , vt ) + 1/2 ) - yt) + *(wr> %)] 

-\ [-*(yJ \vt) + *(2/ 2 + ",2/D ~,yt) + *(y 2 ~ ;vD] (b.i) 



-log 



* Qi r(Q 2 -f.g(y+ + -^-y 2 + --V)) 



vT az vi" vjv2 ViVi 



lQ2T{Ql + y {yt + ^ _ yt _ ^ i i _ _^ j y y - yt 



where 



, T ,/ \ . I dw 2 i , r(i + q (xi + Xl 1 - w 2 -wz 1 )) 

w[x-\,X2) = i ® log t 1 — , (B.2 

V U 2 > f Cl 2ni w 2 - x 2 & T(l - q( Xl + Xi 1 - w 2 - w^ 1 )) ' V ' 

and for §(x\ ) X2) we just note that it starts in any kinematics at least with g 2 . We calculate the 
0(1) expansion of the phase (|B.2|) . Using the property ^>{xx,X2) = ^(xi,0) — ^f(xi, x^ 1 ), being 
valid if \x%\ 7^ 1, and that for \x 2 \ > 1 it starts at g 2 , it is easy to see that we need to calculate 
^(x\,x 2 ) = ^(.TijO) for |jca| < 1, i-e. for x 2 . Since we are interested in the derivative of the expanded 



15 The representations Tl a and ilj, need not be irreducible representations of G. Indeed, in the AdS/CFT case, they 
are sums of irreducible representations, as in I IA.31 . 



3G 



functions with respect to the first argument d\, we need to expand — i(^(j/^~,0) + 1 $'(y 1 , 0)) only. 
Rescaling the integration variable u>2 by g and evaluating the leading residue for small g, we obtain 



^(y+,0)=ilog 



r(i + m 1 + iQ 1 )) 



9{yi ,0) = i log 



r(i + i(pi-iQi)) 



r(i- jfo+iQO) W1 ' °r(i-§(p 1 - l g 1 )) 

The logarithmic derivative of the whole dressing phase is then 



--d Pl io g s*^(p 1 ,p 2 ) = 7 i- 

717 Z7T 



■ V»(i - |(pi + + ^(i + g " ft) + <5i + Q2)) + c.c 



(B.3) 



(B.4) 



where c.c. denotes complex conjugate, and we used that i/)(—i(jp — iQ)) + c.c = ip(l — |(p + iQ)) + c.c 
for integer Q. 

B.2 Performing the integrals by residues 

We demonstrate here how we performed the integrals by evaluating (|5.2ip . In view of the result 
(|B.4[) , we start by evaluating the term with ip(l — \{ftx + iQi)) + c.c. Its contribution factorizes for 
the indices 1, 2 into a product of two factors. The more complicated factor is 



Qi 



dpi 



_ x 2n (pi + Q?U 



V>(i - l(Pi + iQi)) + V(i + \{h - *Qi)) 



(B.5) 



Let us analyze the pole structure of the integrand. Additionally to the two "kinematical" poles at 
p = ±iQ, the polygamma function has poles for ip(—n) if n > 0. These poles are located on the lower 
half plane (LHP) for the first and on the upper half plane (UHP) for second polygamma function: 



-(Q 1 + 2q=ipi) = -n — > pi = Ti(2{n + l) + Q 1 ). (B.6) 

We now use the trick in |62| of exploiting the reality of the integrand to rewrite the integral as 

dpi 



23?e^Q 1 



1 

2n (Pi + Q\) L 



(B.7) 



and close the contour on the UHP. In so doing, we have to pick up the residue at pi = iQi only: 

)L-1 



Q.y- p2 V>(i-|(pi + *Qi)) 



Qi 



J2q 

Q 2 =l 



EE 

Qi j=0 



1 



1-2L+2 



(£-1)! (Pi + iQi) L 

L-l 

( L+ -'~ i ) 

(L-l-j)l 

We now note that the ^-dependent terms give 
dp2 



(-Qtf-L-itfL-i-l) 



(Qi + 1) 



_ 2tt (Pl + Ql) L 



E 



2L - 2 
L - 1 



Hence, the factorizing contribution to , which we denote by A^ , is given by 



Q3-2L 2 1 - 2L = 2 1 ~ 2L 

(1) 



2L-2 
L - 1 



,2L-3 



(B.9) 



L — 



1 jC2L-3EE 

Qi 3=0 



L- 1 St.- - \ 2- 2L + 2 



L + j-1 



(i-i-i)! 



(Ql + 1) 



(B.10) 
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(2) 

Let us concentrate now on the nonfactorizing contributions, which we denote by A s . Using again 
the reality trick, we can write 

E <#a / f / g M + oi)^ + w * (■ + «* ~ «)) ■ 

(B.H) 

and close the p 2 integration contour on the UHP. By picking up the only residue at p 2 = iQ 2 y the 
result is 

q-^2 d V_ V>(i + j(Qi+Q2-i(P2-Pi))) | 
2l Z^ 2 (L-l)! (p 2 + zQ 2 ) i ^ =4Q2 

V2 

= E E( L + r 1 ) 77 ? ^5v(-^) 9 " J; "^ (L " i " 1, (0» + 1 + i(0i + <Pi))-(B.i2) 

Q 2 i=o V J J \^ L J)- * 
The next integral we close on the lower half plane and pick up the residue at — iQ%\ 



L-l 



Qi,Q2h,te=0 v 



-2L+2 



j . • -, \ o-2L+l 

L + J - L ^ rxT(-Q 2 ) 2 - L - J '^ (2i - J1 - J ' 2 - 2) (Qi + Q2 + 1) • (B.13) 



/ (L - l - j 2 ) 

Adding the two terms As = + A^ gives the result we presented in (|5.23[) 
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